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Seismic  Hazards  Studies  for 
Minuteman  Missile  Wings 


1.  INTRODUCTION 

On  28  March  1975,  an  earthquake  of  magnitude  6.  1  occurred  in  the 
Pocatello  Valley  of  southeastern  Idaho.  Although  located  nearly  700  km  east  of 
the  earthquake  epicenter  and  well  outside  of  the  limits  of  the  felt  area  for  this 
event,  the  90th  Strategic  Missile  Wing,  located  at  F.  E.  Warren  AFB,  experienced 
significant  operational  status  changes  as  a  result  of  the  seismic  motions  generated 
by  this  event.  This  situation  indicated  the  need  for  an  examination  of  the  seismic 
hazard  for  each  of  the  Minuteman  Missile  Wings.  This  report  discusses  the 
seismic  hazard  evaluations  conducted  for  the  six  Minuteman  Wings  located  at 
Malstrom,  Ellsworth,  Minot,  Whiteman,  F.  E.  Warren  and  Grand  Forks  Air  Force 
Bases  (Figure  1).  For  each  of  these  facilities,  annual  risk  curves  for  intensity 
and  peak  acceleration,  velocity,  and  displacement  are  presented.  The  uncertain¬ 
ties  which  effect  these  evaluations  are  discussed  and  horizontal  design  response 
spectra  based  on  these  curves  are  given.  Finally,  a  brief  discussion  of  the  local 
seismic  history  for  each  wing  is  presented. 


(Received  for  publication  9  September  1980) 
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Figure  1.  Map  of  the  North  Central  United  States  Showing 
the  Location  of  the  Six  Minuteman  Wings  Studied  in  This 
Report 

2.  SEISMIC  HAZARD  ANALYSIS 
2.1  Seismic  Hazard  Estimation 

Seismic  hazard  analysis  can  be  approached  in  two  wavs;  using  either  deter¬ 
ministic  or  probabilistic  methods.  Deterministic  analysis  relies  on  the  detailed 
knowledge  of  earthquake  faults  in  the  vicinity  of  site  of  interest.  For  each  fault, 
a  maximum  credible  earthquake  is  estimated  and  the  resulting  ground  motions  at 
the  site  of  interest,  using  conservative  assumptions,  are  estimated  for  each  of 
those  events.  This  method  provides  an  upper  bound  estimate  for  the  ground  motion 
at  the  site  of  interest  with  no  consideration  of  the  frequency  of  occurrence.  The 
probabilistic  approach  uses  the  historical  seismicity  in  the  vicinity  of  the  sue  to 
generate  statistical  estimates  of  future  activity,  probabilistic  estimates  of  ground 
motion  at  the  site  can  then  be  calculated  from  this  information.  The  deter:-  •  r t i . 
approach  requires  geologic  field  studies  which  are  typically  not  uv.el  tide  p;,. 
regions  examined  in  this  report.  For  this  reason  only  the  probabilistic  him,.'  lure 
was  applied  in  the  (-valuation  of  seismic  hazards  for  the  missile  u  ings. 

2.2  Seisinicit*  Studies 

The  first  step  in  the  determination  of  probabilistic  hazard  estimates  is  tin- 
development  of  a  statistical  model  of  the  seismic  activity  in  a  region  around  the 
site  of  interest.  Using  historical  earthquake  data,  areas  of  coneent rated  s<  isiec. 
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activity  or  common  tectonic  setting,  known  as  seismic  source  regions,  are  iden¬ 
tified  and  delineated.  All  activity  outside  of  these  regions  is  treated  as  a  unit 
called  the  background  activity.  As  there  is  no  geophysical  basis  for  limiting  the 
location  of  the  background  activity,  these  events  are  assumed  to  occur  with  equal 
probability  throughout  the  area  with  the  exception  of  the  specified  source  regions. 

It  should  be  noted  that  the  lack  of  sufficient  geological  studies,  imprecise 
epicenter  location,  and  incomplete  earthquake  reporting  make  the  boundaries  of 
all  source  regions  somewhat  arbitrary.  For  source  regions  at  large  distances 
front  the  site  of  interest,  the  effects  of  this  uncertainty  are  negligible.  As  the 
site-to-source  distance  becomes  small,  however,  this  problem  can  introduce  sig¬ 
nificant  variability  into  the  seismic  hazard  estimation. 

For  each  of  the  defined  source  zones  and  the  background,  estimates  of  the  level 
of  seismic  activity  were  made  by  fitting  the  standard  recurrence  function 

log  (N)  =  A  -  bM  ( 1 ) 

where 

N  =  the  number  of  events  per  year  of  magnitude  M  or  greater, 

and 

A  and  b  =  the  regression  parameters. 

Prior  to  evaluation  of  the  recurrence  curve  parameters,  the  earthquake  catalogue 

for  each  source  region  was  evaluated  for  temporal  stability  using  a  method  develape 
2 

bv  Stepp.  Where  required  to  obtain  stability,  appropriate  procedures  were  applied 
to  the  earthquake  catalogues.  These  involved  the  editing  of  the  data  o,'»r  selected 
time  windows  for  different  magnitude  ranges  to  produce  stable  estimates  of  re¬ 
currence  periods  within  each  magnitude  range. 

The  uncertainties  associated  with  the  seismicity  study  are  impossible  to  deter¬ 
mine.  While  the  standard  error  of  fitting  for  the  parameters  of  Fq.  (1)  to  the  data 
can  be  determined  easily,  this  is  only  one  source  of  potential  error.  Studies  of 
the  seismicity  of  regions  with  long  historic  earthquake  records,  such  as  China  and 


1.  Kichter,  C.  ( 1 0581  momentary  Seismology,  W,  II.  Freeman  and  t  o., 

San  Francisco,  C  A . 

2.  Stepp,  J.  (1972)  Analysis  of  completeness  of  the  earthquake  sample  in  the  Puget 

Sound  area  and  its  effect  on  statistical  estimates  of  earthquake  hazards. 
Proceedings  of  the  International  Conference  on  Microzonation  for  Safer 
Construction  Research  and  Application,  Seattle,  U'A. 


the  Middle  East,  indicate  that  the  relatively  short  historic  record  in  the  L'nited 

States  is  probably  not  adequate  to  provide  a  complete  description  of  the  seismic 
3  4 

activity.  ’  These  studies  have  shown  that  the  level  and  location  of  seismic 
activity  can  vary  over  periods  longer  than  that  of  thc>  historic  record  in  North 
America.  Apparently  inactive  faults  can  return  to  activity  and  active  regions  be¬ 
come  quiescent.  One  prominent  example  of  the  former  case  is  the  San  Fernando 
Earthquake  of  197  1  which  occurred  on  a  fault  which  appeared  to  be  inactive.  The 
uncertainty  associated  with  this  condition  is  indeterminant  but  could  be  expected 
to  affect  both  the  source  regions  bounds  and  recurrence  curves  used  in  the  risk 
ana  lysis. 

A  maximum  magnitude  earthquake  for  each  source  region  was  also  evaluated 
as  part  of  the  seismicity  studies.  The  magnitude  of  this  event  was  typically  deter¬ 
mined  by  adding  0.  5  magnitude  units  to  the  largest  event  recorded  within  the  source 
area  and  then  rounding  up  to  the  next  higher  integral  or  half  magnitude  value.  This 
value  was  used  as  the  upper  limit  earthquake  in  the  probability  calculations.  While 
the  method  is  arbitrary,  it  is  not  usually  critical  to  the  final  risk  estimation  because 
the  inter-occurrence  periods  of  these  eventsare  typically  much  longer  than  the 
lowest  level  of  risk  examined,  and  contribute  only  slightly  to  the  estimation  of 
ground  motion  at  the  site  of  interest. 

For  each  base  examined  in  this  report,  the  seismic  activity  within  1000  km  of 
the  site  was  used  in  the  seismicity  study.  The  source  of  the  historical  earthquake 
data  was  the  National  Oceanic  and  Atmospheric  Administration  Earthquake  Data 
File.  This  file  covers  earthquake  activity  in  the  United  States  from  1683  to  1978 
but  with  varying  degrees  of  completeness.  Prior  to  1899  this  listing  is  restricted 
to  earthquakes  of  epicentral  intensity  V  or  greater.  Forearthquakes  of  magnitudes 
5.0  rnb  or  greater,  completeness  is  reached  in  the  early  20th  century  for  the  areas 
studied  in  this  report.  Apparent  stability  exists  for  the  range  of  4.  0  to  5.  0  m^  in 
the  early  1960  s  with  the  advent  of  the  World  Wide  Standard  Seismograph  Network 
(WWSSN).  In  general,  seismic  activity  below  magnitude  4.  0  m^  is  reported  only 
erratically.  The  accuracy  of  epicenter  determinations  is  several  tenths  of  a  degree 
in  most  cases. 


3.  McGuire,  K.  (1979)  Adequacy  of  simple  probability  models  for  calculating  felt 

shaking  hazard  vising  the  Chinese  earthquake  catalog,  Hull.  Seism.  Soc.  Am. 
6^:877-892.  . 

4.  Quittmeyer,  1{.  ,  and  Jacob,  K.  (1979)  Historical  and  modern  seismicity  of 

south-central  Asia,  Bull.  Seism.  Soc.  Am.  69:773-823. 
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5.  Meyers,  H.  ,  and  von  Hake,  C.  ( 19  70)  Earthquake  Data  File  Summary, 

National  Geophysical  and  Solar-Terrestrial  Data  Center  Report  KG  HD-5. 


2.3  Ground  Motion  Attenuation 


To  relate  the  temporal  and  spatial  properties  of  seismic  activity  to  ground 
motion  at  the  site  of  interest,  it  is  necessary  to  have  a  means  of  predicting  ground 
motion  at  the  site  as  a  function  of  arthquake  magnitude  and  epicentral  distance. 
Numerous  empirical  ground  motion  attenuation  functions  have  been  developed  for 
this  purpose.  A  typical  equation  has  the  functional  form: 

In  =  a  i  +  a^M  -  a^ln  (K  +  R  )  (2) 

\s  he  re 

g  .  =  the  site  ground  motion, 

M  =  causative  event  magnitude, 

R  =  epicentral  distance, 

H  -  a  distance  factor  commonly  set  to  zero  or  25  km, 

and 

-ij.  a a =  empirically  derived  constants. 

The  groum.  motion  descriptor,  g^,  is  usually  peak  site  acceleration,  velocity,  or 
displacement;  though  other  parameters  have  been  used. 

Virtually  all  of  the  available  attenuation  functions  have  one  significant  deficiency. 
They  have  been  derived  on  the  basis  of  data  primarily  from  Southern  California 
earthquakes.  The  ground  motion  attenuation  in  this  region  is  not  thought  to  be 
typical  of  the  rest  of  the  United  States.  For  this  reason,  these  equations  can  not 
be  applied  blindly  outside  of  California. 

To  overcome  this  problem,  two  different  paths  were  followed  to  obtain  ground 
motion  estimates  at  the  various  wings.  The  first  of  these  was  to  use  Modified 
Mercalli  intensity  (Table  1)  as  the  ground  motion  descriptor.  As  attenuation  curves 
for  intensity  can  be  derived  on  the  basis  of  historical,  non-instrumental  earthquake 
records,  they  are  generally  available  for  all  of  the  United  States.  1'he  difficulty 
with  this  approach  is  the  subsequent  conversion  from  intensity  to  ground  motion  des¬ 
criptors  of  engineering  significance.  As  the  intensity  scale  was  based  on  subjective 
and  coarse  gradations,  wide  ranges  of  acceleration,  velocity,  and  displacement 

g 

have  been  observed  for  each  intensity  level  (Figure  2).  The  conversion  from  site 

6.  Trifunac,  M. ,  and  Brady,  A.  ( 1975)  On  the  correlation  of  seismic  intensity 

scales  with  peaks  of  recorded  strong  ground  motion.  Bull.  Seism.  Soc.  Am. 
65:139-  162. 
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intensity  to  peak  acceleration,  velocity,  and  displacement  were  made  using  the 
relationships  derived  by  Trifunac  and  Hrady.  ^  Specifically: 

log  =  -0.  014  +  0.  301  ,  (3) 

b  H  mm 

log  V„  =  -0.  63  +  0.  251  ,  (4) 

&  11  mm 

and 

log  D..  =  -1.  13  +  0.  241  ,  (5) 

6  H  mm 

where 

A..,  V...  Du  =  peak  horizontal  acceleration,  velocity,  and  displacement 
M  rl  ri  ‘ 

respectively,  and 


1  =  the  site  Modified  Mercalli  intensity, 

mm  J 

These  relations  were  explicitly  derived  over  the  range  of  intensities  from  IV  to  X 
but  have  been  assumed  to  apply  over  the  full  twelve-point  range  of  the  scale  for  the 
purposes  of  this  report.  It  should  also  be  noted  that  the  use  of  other  empirical 
relations  can  produce  different  ground  motion  levels. 

The  second  approach  that  was  used  to  develop  the  seismic  risk  estimates  was 
to  attempt  to  derive  regionally  modified  peak  ground  motion  attenuation  functions 
for  the  areas  being  studied.  Only  for  acceleration  was  this  procedure  successful. 

A  description  of  the  method  used  for  regional  modification  of  peak  acceleration 
attenuation  functions  is  given  in  Appendix  A.  The  results  of  this  analysis  indicate 
that  only  in  the  case  of  Wing  IV,  Whiteman  Al'H,  Missouri,  is  the  regional  attenua¬ 
tion  so  different  from  that  in  California  as  to  require  an  attenuation  function  other 
than  an  empirical,  California  data  based  equation.  Since  the  completion  of  this 
study,  several  useful  suggestions  to  improve  the  modification  scheme  of  Appendix  A 
have  been  received  (Xuttli,  personal  communication,  lhtiO).  Tile  principal  effect 
of  these  changes  would  be  to  alter  the  near-field  acceleration  estimates  for  the 
Central  United  States  attenuation  function.  The  hazard  estimates  at  all  of  the  sites 
except  Whiteman  AT'B  are  not  expected  to  be  affected.  Preliminary  results  suggest 
a  reduction  of  the  estimated  seismic  risk  at  Whiteman  Al'H,  however,  the  degree 
of  change  is  uncertain.  At  the  present  time,  the  modifications  are  being  incorpora¬ 
ted  into  the  technique  and  a  new  Central  United  States  attenuation  function  evaluated. 
After  this  has  been  completed  the  seismic  hazard  at  Whiteman  Al'H  will  be  re¬ 
evaluated  and,  if  significant  changes  result,  the  new  risk  estimate  will  be  pre¬ 
sented  at  a  later  date. 
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Table  1.  Modified  Mercalli  Intensity  Scale  of  1931 


1.  Not  felt  except  by  a  very  few  under 
specially- variable  circumstances. 
(I  Kosse-Forel  scale.  ) 

U.  Kelt  only  by  a  few  persons  at  rest, 

especially  on  upper  floors  of 
buildings.  Delicately  suspended 
objects  may  swing.  <1  to  II 
Rossi- Korel  scale. ) 

III.  Kelt  quite  noticeably  indoors, 

especially  on  upper  floors  of 
buildings,  but  many  people  do  not 
recognize  it  as  an  earthquake. 
Standing  motorcars  may  rock 
slightly.  Duration  estimates. 

(Ill  Rossi- Korel  scale.  ) 

IV.  During  the  day  felt  indoors  by  many, 

outdoors  by  few.  At  night  some 
awakened.  Dishes,  windows, 
doors  disturbed;  walls  make 
creaking  sound.  Sensation  like 
heavy  truck  st l  iking  building. 
Standing  motorcars  rocked 
noticeably.  (IV  to  V  Rossi- Korel 
scale. ) 

V.  Kelt  by  nearly  everyone,  many 

awakened.  Some  dishes,  windows, 
etc.  ,  broken;  a  few  instances  of 
cracked  plaster;  unstable  objects 
overturned.  Disturbance  of 
trees,  poles,  and  other  tall 
objects  sometimes  noticed. 
Pend-lum  clocks  may  stop. 

(V  to  VI  Rossi-Korel  scale.  ) 

VI.  Kelt  by  all,  many  frightened  and 
run  outdoors.  Some  heavy 
furniture  moved;  a  few  instances 
of  fallen  plaster  or  damaged 
chimneys.  Damage  slight. 

(VI  to  VII  Rossi-Korel  scale.  ) 

VII.  Kverybody  runs  outdoors.  Damage 
negligible  in  buildings  of  good 
design  and  construction;  slight 
to  moderate  in  well-built  ordina ry 
structures;  considerable  in 
poorly  built  or  badly  designed 
structures;  some  chimneys 
broken.  Noticed  by  persons 
driving  motorcars.  (VIII 
Rossi-Korel  scale.  ) 


V 111.  Damage  slight  in  specially 

designated  structures;  con¬ 
siderable  in  ordinary  sub¬ 
stantial  buildings  with  partial 
collapse;  great  in  poorly  built 
structures.  Panel  walls  thrown 
out  of  frame  structures.  Fall 
of  chimneys,  factory  stacks, 
columns,  monuments,  walls. 
Heavy  furniture  overturned. 
Sand  and  mud  ejected  in  small 
amounts.  Changes  in  well 
water.  Persons  driving 
motorcars  disturbed.  (VIII-+ 
to  IX-  Rossi-Korel  scale.  ) 

IX.  Damage  considerable  in  specially 
designed  structures;  well- 
designed  frame  structures 
thrown  out  of  plumb;  great  in 
substantial  buildings,  with 
partial  collapse.  Buildings 
shifted  off  foundations. 

Ground  cracked  conspicuously. 
Underground  pipes  broken. 

(IX+  Rossi-Korel  scale.  ) 

X.  Some  well-built  wooden  struc¬ 
tures  destroyed;  most 
masonry  and  frame  structures 
destroyed  with  foundations, 
ground  badly  cracked.  Rails 
bent.  Landslides  consider¬ 
able  from  riverbanks  and 
steep  slopes.  Shifted  sand 
and  mud.  Water  splashed 
(slopped)  over  banks.  (X 
Rossi-Korel  scale.  ) 

XL  Few,  if  any,  (masonry) 

structures  remain  standing. 
Bridges  destroyed.  Broad 
fissures  in  ground.  Under¬ 
ground  pipelines  completely 
out  of  service.  Karth  slumps 
and  land  slips  in  soft  ground. 
Rails  bent  greatly. 

XII.  Damage  total.  Waves  seen  on 
ground  surfaces.  Lines  of 
sight  and  level  distorted. 
Objects  thrown  upward  into 


15 


r 


10  t- 


10  \r 


*0  fc" 


Accelerotior 

(cm/sec2) 


Velocity 

(cm/sec) 


fl 


Displacement 
(cm  ) 


Figure  2.  Mean  and  One 
Standard  Deviation  Error 
Bars  Peak  Ground  Motions 
Plotted  Against  Modified 
Mercalli  Intensity  at  the 
Observation  Site.  Accelera¬ 
tion,  open  circles;  velocity, 
open  squares;  displacement, 
solid  circles 
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If  the  analysis  of  Appendix  A  is  correct,  however,  it  would  appear  that  most 
empirically  derived  functions,  being  based  predominantly  on  data  from  less  than 
100  km,  over-estimate  far-field  accelerations.  In  Figure  3,  the  predicted  peak 
accelerations  for  a  magnitude  G.  0  earthquake  in  California  have  been  plotted  using 
the  attenuation  functions  derived  in  Appendix  A  and  seven  empirical  functions  listed 
by  McGuire.  1  The  diserepeney  at  distances  greater  than  about  100  km  is  apparent. 
This  problem  is  significant  in  th.it  for  most  sites  studied  in  this  report,  the  seismic 
source  regions  tend  to  be  at  distances  greater  than  100  km.  In  turn,  one  must  also 
be  skeptical  of  the  velocity  and  displacement  relations  for  the  same  reason.  The 
set  of  velocity  and  displacement  attenuation  functions  which  predict  the  lowest  mean 
values  at  large  opiecntr.nl  distances  were  those  developed  by  Orphal  and  l.ahoud.  ^ 
these  equations  were  used  for  risk  estimation  at  all  wings  except  Wing  IV.  The 
constants  of  all  ground  motion  attenuation  functions  used  in  this  report  are  given 
in  Table  2. 


7.  McGuire,  11.  (197(i)  iOKTHAN  Computer  Frog  ram  for  Seismic  disk  Analysis, 

F.S.  Geol.  Sure.  Open-File  Kepoi-t  TG-GT. 

8.  Orphal,  D.  ,  and  I  .a  hour],  .1.  ( 1 !  >  <  4 )  Prediction  of  peak  ground  motion  from 

earthquakes,  Hull.  Seism.  Sec.  Am.  fid:  1  aG3-  1  fi 7 4. 
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Figure  3.  Plot  of  Peak  Ground  Acceleration  vs  Distances  for  a  Magnitude 
6.  0  Ml  Harthquake  in  Southern  California  as  Predicted  by  Seven  Kmpiricul 
Curves  (dashed  lines)  and  the  Derived  Curved  (solid  line)  Based  on  the 
Analysis  of  Appendix  A 


Table  2.  Ground  Motion  Attenuation  Functions 


Cl  round  Motion 

Application 

Source 

ln  Ks 

1 

2 

£ 

M 

In  (It  * 

K  ) 
o 

al 

lV3 

K 

o 

a, 

lo 

Intensity  ^ 

1,  II,  111, 

V,  VI 

Howell  aruj 
Schultz ' 

2.  2!) 

1.  -42 

I.  36 

0.  0 

0.  68 

Intensity  ^ 

IV 

Howell  and 
Schultz" 

I.  20 

2.  23 

1. 36 

0.  0 

0.  68 

Acceleration 

I.  II,  III, 

V,  VI 

Appendix  A 

8.  76 

0.  8  12 

2.  32, ' 

25.  0 

(0.  7  07  B 

Ac  coloration 

IV 

Appendix  A 

■4.21 

1.  -40 

1. 57 

25.  0 

(0.  707  1  :■ 

Velocity 

I.  II.  Ill, 

V.  VI 

Orphal  and 
I.ahoud® 

-0.  32 

1.  20 

1.  3-4 

0.  0 

(0.  7  07  )* 

1  )is  placement 

I,  II,  III, 

V,  VI 

Orphal  and 
Lahoud® 

-3.  06 

1.312 

I.  18 

0.  0 

(0.  7  601* 

^Conversion  from  I  lo  M  using  relationships  bv  Brazen** 


Assuined  value. 
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2.4  Seismic  Risk  Estimation 


9 

Using  a  method  proposed  by  Cornell'  and  implemented  in  a  FORTRAN  com - 
7 

outer  program  by  McGuire  the  temporal  and  spatial  distribution  of  seismic 
activity  and  the  ground  motion  attenuation  fun.  tions  ean  be  t-omirn  --1  int,,  sing!,- 
statement  of  the  probability  of  attaining  a  given  level  of  ground  motion  at  the  site 
of  interest.  In  this  procedure,  the  probability  that  the  ground  motion  level  will 
reac  h  or  exceed  a  specified  level,  m^,  is  defined  as  the  integral  of  the  product  of 
the  independent  probability  density  functions  for  magnitude,  fy,  distance,  t‘  and 
distance,  r.  This  can  be  stated  as  an  equation: 


P[M  *  m  |  --  j  J  P[M  2  m  |  s  and  r|  f  (s)  f  (r)dsdr 
fs  s  H  &  o  r\ 


(6) 


whore 


P[M 
1  g 
distance. 


m 


g 


and  r|  is  the  condition  probability  given  event  magnitude  and 


The  conditional  probability  of  Eq.  (6)  is  a  function  of  the  ground  motion  attenua 
tion  equation  and  its  standard  deviation.  The  function  f^is)  is  derived  from  each 
source  region  recurrence  curve  while  f^(r)  incorporates  the  spatial  relationship 
between  source  region  and  site  of  interest.  Evaluation  of  the  integral  debts  the 
probability  of  one  event  from  the  specified  source  region  reaching  or  -ceding  n 

i 

Bv  multiplying  this  value  bv  the  expected  number  of  .-vents  in  -be  region  m  l 
accumulating  nv  -r  all  source  areas,  the  total  expected  number  of  events  .-.-ting 
the  conditions,  mj.  is  obtained.  Assuming  earthquake  . . .  ,  urr>n<  -  is  a 

Poisson  process,  the  annual  risk  is  given  hv 

K[  M  ’  m  | 

RIMr  '  "\,l  '  -  e  U  (7) 

Annual  risks  -  an  be  .-.inverted  to  risks  for  anv  given  lifetime  bv  the  equation; 

»N  1  -  «a»N  (») 


where 


N  is  the  lifetime  in  years 
R  ^  is  the  annual  risk, 


and 


R  is  the  N'-vear  lifetime  risk. 

-\ 

One  final  term  used  in  describing  risk  is  the  return  period  which  is  the  reciprocal 
of  the  risk  and  describes  the  number  of  lifetime  periods,  on  the  average,  between 
recurrence  of  the  associated  level  of  ground  motion. 

2.5  Response  Spectra 

For  engineering  purposes,  the  spectral  characteristics  of  ground  motion  arc 
tvpicallv  presented  in  the  form  of  response  spectra.  These  spectra  represent  the 
maximum  response  of  a  simple,  viscous -dam ped  harmonic  oscillator  over  a  range 
of  natural  periods  for  a  specified  percentage  of  critical  damping.  Methods  have 
been  developed  to  estimate  upper  limit  response  spectra  given  the  level  of  ground 
motion  at  the  site  of  interest,  known  as  design  response  spectra.  The  AKC 
amplifications  factors,  given  in  Table  3,  are  used  to  modify  the  predicted  ground 
motion  levels  at  a  site  to  obtain  the  response  spectra  levels  at  the  specified 
frequencies.  These  factors  give  the  mean  plus  one  standard  deviation  response 
spectra. 


Table  3.  Horizontal  Design  Response  Spectra  Amplification  Factors 


' 

Acceleration 

Displacement  (cm) 

Critical  Damping 

(%> 

33  Hz 

0  Hz 

2.  5  Hz 

0.  25  Hz 

0.  5 

1.  0 

4.  06 

5.  05 

3.  20 

2.  0 

1. 0 

3. 54 

4.25 

2.  50 

5.  0 

1.0 

2.  61 

3.  13 

2.  05 

7.0 

1.  0 

2.  21 

2.72 

1. 88 

10.  0 

1.  0 

1.  01) 

2.  28 

1.70 

The  levels  of  critical  damping  correspond  to  soil  conditions  at  the  site  with 
lower  values  associated  with  harder  ruck.  Most  foundation  conditions  lie  between 
2.0  and  7.0  percent  of  critical  damping.  The  lowest  damping  is  always  the  higher 
amplitude  spectrum.  In  all  response  spectra  plots  m  this  report  the  complete 


10.  Hays.  \V.  ,  Algcrmissen.  S.  ,  Kstinesi.  A,.  I’erkms,  I).,  and  Rinehart.  W. 
( 1 1 '7 a )  f  i ui del i ivs  for  lievclopmg  Design  Earthquake  Response  Sped  ra, 

I ' .  S.  ( 1  ciil .  iTii rv .  TocTm i  cal  Report  \l  -  111. 
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range  of  damping  factors  given  in  Table  3  arc  plotted.  Thus,  the  upper  plot  is 
always  0.5  percent  and  the  lower  10  percent  of  critical  damping. 

It  should  be  noted  that  when  the  design  spectra  is  calculated  on  the  basis  of 
independently  evaluated  acceleration  and  displacement  levels  corresponding  to  a 
specified  risk  level,  they  are  not  necessarily  a  spectral  limit  for  any  one  earth¬ 
quake.  This  is  because  the  ground  motions  on  which  they  are  based  may  be 
generated  by  different  earthquakes-.  *  *  In  this  case,  it  is  more  accurate  to  view 
the  spectra  as  upper  limits  for  separate  frequency  bands  than  as  the  upper  limit 
spectra  for  any  one  earthquake.  The  spectra,  in  this  case,  will  be  designated  as 
composite  response  spectra. 

3.  SEISMIC  HAZARD  EVALUATIONS 
3. 1  Malstrom  AFB  —  Wing  I 

Malstrom  AFB  is  located  near  one  of  the  most  seismic-ally  active  regions  of 
the  continental  United  States,  an  area  known  as  the  Intermountain  Seismic  Belt. 

This  zone  runs  along  the  northern  Rocky  Mountains  in  western  Montana  through 
Yellowstone  National  Park  and  along  the  Wasatch  Front  in  Utah.  The  location  of 
Malstrom  AFB  relative  to  those  active  regions,  as  shown  in  Figure  4,  suggests  a 
high  level  of  seismic-  hazard  at  this  facility. 

For  the  purposes  of  this  studv,  the  Intormountain  Seismic:  Belt  has  been 
divided  into  three  subregions,  Montana,  Yellowstone  and  Utah,  to  compensate  for 
the  variations  of  seismic  activity  tin  oughout  the  zone.  In  addition,  several  other 
source  regions  were  defined  bv  analysis  of  epicentral  clustering  and  tectonic  setting 
of  all  earthquakes  within  1000  km  of  Malstrom  AFB.  These  source-  regions  are 
depicted  in  Figure  4.  The  area,  recurrence  curve  constants  and  maximum  magni¬ 
tude  earthquakes  for  cacti  source  region  are  given  in  Table  4  along  with  the  back¬ 
ground  seismicity  parameters.  The  background  recurrence  curve-  is  assumed  to 
hold  in  all  areas  not  explicitly  in  a  seismic-  source  zone'.  Thus,  while-  the-  prob¬ 
ability  of  an  earthquake  occurring  outside-  of  a  defined  source  is  low,  it  is  assumed 
that  activity  can  occur  anywhere  within  1000  km  of  the  site-  o;  interest  including 
areas  with  no  known  historical  reports  of  activity.  The-  site-  used  to  i-e-pre-se-nt 
Malstrom  AFB  in  the  risk  calculations  was  47.  5  ;  X  and  1  I  1.  0°w. 

11.  Battis,  .1.  ( 1078a)  (ie'ophvsical  Studies  for  Missile  Basing  Se  ismic  Risk  Studies 
in  the-  Wi-stern  lYiit  e-d  .States,  Texas  Instruim-nts  Inc.,  ?>ci'-iit  iTTc  Report  No.  Li 
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Figure  4.  Seismic  Source  Regions  Used  to  evaluate  the 
Seismic  Hazard  at  Walstrom  AFB  (indicated  by  solid 
circle) 


Table  4.  Malstrom  AFB  Source  Region  Parameters 


Area 

4  9 

MAX 

Source 

(10  km  “1 

A 

B 

Mc 

Wyoming  -  Colorado  - 
South  Dakota 

43.05 

mm 

0.  592 

7.  0 

Nevada 

16.  66 

m&M 

1. 031 

7.0 

Washington  -  Oregon 

26.  66 

2.  806 

0.  597 

7.0 

Central  Idaho 

4.  41 

3.  177 

0.  792 

7.  0 

V  ellowstone 

6.  1 1 

3.  652 

0.765 

7.0 

Montana 

7.70 

3.  101 

0.  682 

7. 0 

Utah 

20.  05 

3.  23  6 

0.  690 

7. 0 

Background 

2.  636f 

0.  7  48 

6.  0 

log  ( \  /  v  /  1 0  km  ) 


The  estimated  annual  risk  curve  lor  Modified  Mercalli  intensity  is  shown  m 
Figure  ">  and  the  peak  acceleration,  velocity,  and  displacement  curves  are  displace 
in  Figure  6.  The  ground  motion  levels  at  specified  annual  risks  are  also  tabulated 
in  Table  a.  In  1  aid  e  fi  the  ground  motion  levels  based  on  the  .  envoi's  ion  of  Modi  1 1  e  I 


Mercalli  intensity  are  given.  Considering  the  large  standard  errors  assoejat>  I 


with  the  conversion  of  site  intensity  to  the  peak  ground  motion  levels,  the  values 
in  these  tables  are  believed  to  be  in  good  agreement. 


Figure  5.  Modified  Mercalli  Intensity  Seismic  Hisk  Curve  for 
Malstrom  AFB 


Table  6.  Peak  Ground  Motion  Annual  Risk  Levels  for  Malstrom  AKll 
Based  on  Intensity  Levels 


Annual 

Risk 

Fractional 

Intensity 

Acceleration 
(cm /sec*) 

Velocity 

(cm/sec) 

Displacement 
(cm ) 

0.  9 

1.  8 

3.  5 

0.7 

0.  2 

0.  5 

2.  9 

6.  9 

1.2 

0.  4 

0.  2 

3.  6 

11.  6 

1.9 

0.  5 

0.  1 

4.  0 

15.  3 

2.  3 

0.  7 

0.  05 

4.  4 

19.  5 

2.  9 

0.  8 

0.  02 

4.  8 

26.7 

3.7 

1.  1 

0.  01 

5.  1 

32.  8 

4.  4 

1. 2 

0.  005 

5.  4 

39.  8 

5.  1 

1.  5 

0.  002 

5.7 

49.  6 

6.  2 

1.7 

0.  001 

5.  9 

57.  8 

7.  1 

2.  0 

The  proximity  of  the  Montana  source  region  to  the  site  of  interest,  as  shown 
in  Figure  4,  suggests  that  a  strong  spatial  gradient  of  the  seismic  hazard  should 
be  observed  in  the  vicinity  of  Malstrom  AFB.  This,  in  fact,  is  the  case  with 
sharply  increasing  ground  motions  at  each  level  of  risk  observed  toward  the  south¬ 
west  of  the  facility.  In  Figure  7,  the  (teak  ground  accelerations  for  an  annual  risk 
of  0.002,  or  return  period  of  500  years,  arc  contoured  over  a  one-degree  square 
centered  on  Malstrom  AFB.  For  other  levels  of  risk  the  contours  follow  in  the 
same  trend.  The  same  characteristics  are  observed  for  peak  velocity  and  dis¬ 
placement  although  tiie  gradient  is  less  severe  than  for  acceleration. 

It  should  be  noted  that  the  contour  locations  arc  defined  bv  the  location  of  the 
Montana  source  region  boundaries.  Inside  of  the  boundary,  neglecting  the  effects 
of  other  source  regions,  the  seismic  hazard  is  uniform.  As  the  boundary  is 
approached  what  could  be  termed  an  edge  effect  begins  to  reduce  iho  hazard. 

Moving  awav  from  the  source  the  hazard  decreases  gradually  :o  a  uniform  level 
determined  by  tin1  background  activity.  Moving  the  source  region  boundary  towards 
the  site  of  interest  will  increase  the  seismic  hazard  at  the  site  and,  conversely, 
mov  ing  the  boundary  awav  from  the  site  will  reduce  the  hazard.  'I  lie  sourc  re¬ 
gion  boundaries  arm  somewhat  arbitrary  due  to  a  lack  of  detailed  fault  maps  of  a 
complete  un  li  rstanding  of  the  geologic  setting  associated  with  earthquakes  m  this 
region.  For  sites  near  source  regions,  this  rail  produce  a  decrease  in  tile  level 
of  confidence  associated  with  the  estimated  ground  motion  values. 
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Figure  7.  Contours  of  500-vear 
Return  Period  Accelerations  in 
a  One-degree  Square  Centered 
on  MaLstrom  AFB.  (Solid 
circle);  (units  in  7,V.  g) 


In  Figures  8a,  8b,  and  8c  the  composite  horizontal  design  response  spectra 
for  10-.  100 -,  and  1000 -year  return  period  ground  motions,  based  on  the  values 
given  in  Table  5,  are  displayed. 

In  Figure  0,  95  earthquakes  between  1972  and  1978  reported  in  the  Earthquake 
Data  File  within  200  km  of  Malstrom  AFB  are  plotted.  The  closest  center  of  major 
activity  to  the  site  is  120  km  to  the  southwest  near  Helena,  Montana  although  an 
earthquake  with  no  recorded  magnitude  or  intensity  was  reported  within  75  km  of 
tiie  base.  The  largest  event  reported  in  this  particular  area  was  a  magnitude 
6.25  M  ^  event  which  occurred  on  19  October  1935  and  is  one  of  the  largest  earth¬ 
quakes  reported  in  the  northern  extent  of  the  Intermountain  Seismic  Belt.  This 
earthquake  was  one  of  a  series  of  events  occurring  in  late  1935  of  which  four  events, 
including  the  magnitude  6.  25  M^  earthquake,  had  epieentral  intensities  of  between 
VI  and  V  III.  At  the  90  percent  confidence  level  the  predicted  ground  motions  at 
Malstrom  AFB  lor  this  event  are  21.9  on/soc",  5.  1  em/soc  and  1.5  em/sec.  The 
horizontal  design  response  spectra  for  this  event  is  shown  in  Figure  10.  An  aver¬ 
age  site  intensity  determined  from  these  values.  Modified  Meroalli  Intensity  Y, 

1 

agrees  with  the  intensity  level  observed  m  the  vicinity  of  Great  Falls,  Montana. 


12.  Neumann,  F.  (1937)  United  States  Earthquakes.  1935.  Department  of 
Commerce,  Coast  and  Geodetic  Survey,  Serial  N'o.  6(10. 
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Figure  8a.  Composite 
Horizontal  Design 
Response  Spectra  for 
Malstrom  AFB  Based 
on  10-year  Return  Period 
Ground  Motions 


Figure  8b.  Composite 
Horizontal  Design 
Response  Spectra  for 
Malstrom  AFB  Based 
on  100-year  Return  Period 
Ground  Motions 


Figure  8c.  Composite 
Horizontal  Design 
Response  Spectra  for 
Malstrom  AFB  Based  on 
1000-vear  Return  Period 
Ground  Motions 
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Figure  9.  Reported  Earthquake  Epicenters 
Within  200  km  of  Malstrom  AFB 

At  a  greater  epicentral  distance,  170  km,  an  even  larger  earthquake,  magni¬ 
tude  6.75  m,  was  reported  on  28  June  1925.  The  90  percent  confidence  ground 
motions  at  the  site  of  interest  are  estimated  at  17.2  cm/sec  ,  6.0  cm/sec,  and 
2.0  cm.  The  lower  acceleration  and  higher  velocity  and  displacement  as  compared 
to  the  Helena,  Montana  Earthquake  are  explained  by  the  longer  travel  path  and  the 
resulting  greater  attenuation  of  the  high  frequency  components  of  the  signal.  Again, 
the  predicted  site  intensity  of  this  event  at  Great  Falls,  Montana  is  intensity  V. 

The  horizontal  design  response  spectra  are  shown  in  Figure  11.  The  higher 
amplitudes  at  low  frequency  than  for  the  Helena  event  are  the  significant  feature  of 
this  event. 

3.2  Ellsworth  AFB  —  Wing  II 

In  Figure  12,  the  seismic  source  regions  used  in  Ihe  hazard  evaluation  for 
Ellsworth  AFB  are  shown  and  the  associated  source  parameters  are  given  in 
Table  7.  Several  of  (he  source  regions  have  been  redefined  as  compared  to  Ihose 
used  in  the  previous  section.  This  results  from  the  change  in  the  earthquakes 
included  in  the  seismicity  studies  due  to  the  1000  km  radius  search  limit  from  the 
site  of  interest.  In  addition,  as  a  result  of  the  greater  influence  on  seismic  hazard 


by  near  site  source  regions,  greater  care  was  used  in  defining  the  Colorado  and 
Wyoming  area  source  regions.  The  point  used  in  all  calculations  for  Ellsworth  AFB 
was  44.  2°N  and  103.  0°W. 


Figure  10.  Horizontal 
Design  Response  Spectra 
for  Malstrom  AFB 
Assuming  a  Recurrence 
of  the  1935  Helena,  Montana 
Earthquake  (6.25  M^) 


Figure  11.  Horizontal 
Design  Response  Spectra 
for  Malstrom  AFB 
Assuming  a  Recurrence  of 
the  1925  Earthquake 
(6.  75  M  ) 


Ktlsworth  Al'  li  Snurco  K<\t»ion  Parameters 


S'HlI'ic 

Area 

( 104  km2) 

gglgig 

HH 

m,max 

mgm 

BH| 

Central  Mah" 

•4.41 

3.  177 

0.  792 

7.  0 

Montana 

7.70 

3.  101 

0.  682 

7.  0 

V  el  Low  stone 

6.  1 1 

3.  652 

0.  765 

7.  0 

Colorado 

15.  83 

3.  4  1!) 

0.  700 

7.  0 

Wyoming 

18.  05 

4. 465 

1 .  114 

6.  0 

I’tah 

26.  10 

3.  5!>8 

0.  702 

7 ,  f> 

.Nevada 

7.  58 

1.  662 

0.  5  16 

7 . 0 

New  Mexico  -  Texas 

!>.  15 

2.  57  4 

0.  683 

7.5 

Background 

1.721  4 

0.  199 

6.  0 
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I-igore  12.  Sris  n>  i  r  Source  Hrgions  1'srd  In  Kvnluat 
Seismic  Hazard  at  Kllswnrlh  .m  l  I'.  I).  Warren  Air  I 
Bases  (indicated  liv  dosed  rin  les) 


The  annual  risk  curves  for  Modified  Mercalli  intensity  and  peak  ground  aecelera 
tion,  velocity,  and  displacement  for  Ellsworth  AFB  are  displayed  in  Figures  13 
and  14.  Table  8  gives  the  ground  motions  calculated  at  specific  levels  of  annual 
risk  for  Ellsworth  AFB.  A  second  set  of  peak  acceleration,  velocity,  and  displace¬ 
ment  motion  levels  was  constructed  by  conversion  of  the  estimated  site  intensities. 
These  values  are  given  in  Table  9.  Given  the  limitations  of  calculating  ground 
motions  from  site  intensity,  the  two  sets  of  values  are  considered  in  good  agree¬ 
ment.  This  is  especially  true  in  that  the  site  intensity  values  are  only  mean  values 
at  the  specified  intensities.  For  example,  while  the  0.001  annual  risk  accelera¬ 
tions  appear  to  be  significantly  different  the  value  calculated  directly  is  within  the 
90  percent  confidence  bounds  for  site  intensity  V. 


n  m  n  i 

MODIFIED  MERCALLI  INTENSITY 


Figure  13.  Modified  Mercalli  Intensity  Seismic  Risk  t'urv 
for  Ellsworth  AFB 


Killin'  14.  Peak  G round  Acceleration,  Velocity,  and  Displacement 
Seismic  Risk  Curves  for  KUsworth  AKB 


Table  8.  Peak  Ground  Motion  Annua!  Risk  Bevels  for  KUsworth  AKB 


Annual 

Risk 

Return 

Period 

(years) 

Intensity 

H 

Acceleration 

(cm/sec2) 

Velocity 
(cm  /  sec) 

Displacement 
(cm ) 

0.  9 

1.  11 

I 

1. 

6 

0.  1 

0.  03 

0.  5 

2 

II 

3. 

4 

0.  2 

0.  07 

0.  2 

5 

II -III 

6. 

4 

0.  4 

0.  14 

0.  1 

10 

III 

9. 

3 

0.  6 

0.  21 

0.  05 

20 

III 

13. 

1 

0.  9 

0.2  9 

0.  02 

50 

IV 

20. 

3 

1.  3 

0.  4  5 

0.  01 

100 

IV 

27. 

9 

I.  8 

0.  60 

0.  005 

200 

IV 

38. 

2 

2.  4 

0.  79 

0.  00  2 

500 

V 

57. 

3 

3.  6 

1.11 

0.  00  1 

1000 

V 

77 . 

2 

4.  8 

1.  42 

Table  9.  Peak  Ground  Motion  Annual  Risks  Levels  at  Ellsworth  APB 
Based  on  Intensity  Levels 


Annual 

Risk 

Fractional 

Intensity 

Acceleration 
(cm  'se:  -! 

o.  y 

L  1 

2.  1 

1 

0.  4 

o.  i 

0.  5 

2.  1 

4.  0 

0.  8 

0.  2 

0.  2 

2.  8 

6.  5 

1.  1 

0.  3 

0.  1 

3.  2 

8.  8 

1.  5 

0.  4 

0.  05 

3.  6 

11.4 

1.  8 

0.  5 

0.  02 

4.  1 

15.  8 

2.  4 

0.  7 

0.  01 

4.4 

19.  3 

2.  8 

0.  8 

0.  005 

4.7 

23.  5 

3.  3 

1.  0 

0.  002 

5.  0 

30.  6 

4.  2 

1.  2 

0.  001 

5.  3 

37.  5 

4.  9 

1.  4 

For  this  wing,  the  acceleration  levels  for  annual  risk  are,  to  a  large  degree, 
determined  by  the  background  seismicity  characteristics.  For  peak  velocity  and 
displacement,  which  attenuate  slower  than  acceleration,  t lie  seismic  source  re¬ 
gions  in  Wyoming  and  Colorado  contribute  to  the  risk  in  a  greater  amount.  In 
neither  case,  however,  do  t ho  nearest  sources  play  a  dominant  role  in  establishing 
the  seismic  risk.  In  this  case,  the  spatial  change  in  seismic  hazard  is  expected  to 
be  negligible.  Little  change  in  the  acceleration  values  would  be  expected  within 
50  km  of  the  site  in  any  direction.  Verv  gradual  increases  in  velocity  and  dis¬ 
placement  could  be  expected  towards  the  southwest,  as  the  defined  sources  are 
approached,  but  within  50  km  the  changes  are  not  exported  to  bo  substantial. 

In  Figures  15a,  15b,  and  15c,  the  10-,  100 -,  and  1000 -vear  return  period 
composite  horizontal  design  response  spectra  determined  for  Ellsworth  Al'll  are 
shown.  These  spectra  are  based  on  the  ground  motion  levels  given  in  Table  fi. 

Twelve  moderate  sized  earthquakes  have  been  reported  within  200  km  >1 
Ellsworth  AI  M  between  1895  and  1978.  The  epicenters  of  these  events  are  shown 
in  f  igure  10.  The  .-loses!  event  was  of  magnitude  I.  1  U  |  and  epieentral  intensity 
\  I  loeated  24  km  to  the  west  of  the  base  and  oeeurre.l  on  20  .nine  l''0M.  I  la 
depicted  events  ire  cr rat  i .  all  v  di st  i diluted  m  sp  i.  e  an  i  ;c-o\  i  ie  .  i : ! , ■  1 1 1 fo re  . at  ton 
on  future  potential  epicenters  w  ilhin  tins  region.  I  tr-  r-.axie  ue  :  >  not  t<  I  1  ,  _<m 
tude  in  t his  set  of  event  s  w  as  5.  1  t ■ ,  and  cpi . -  cut  i  d  i  n' "t,  -  r  v  .  !  I ,  !-■•!:  ■  i ■  • ; ■  •  i  t  m.  i 
for  tlie  same  event. 
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Figure  16.  Reported  Earthquake  Epicenters 
Within  200  km  of  Ellsworth  AFB 

3.3  Minot  AFB  and  Grand  Forks  AF'B  —  Wings  HI  and  V| 

As  with  the  previous  studies  a  search  of  the  NOAA  Earthquake  data  file  was 
made  for  all  earthquakes  within  1000  km  of  tiiese  two  facilities.  While  1100 
events  are  located  within  that  radius  of  Minot  AFB,  only  124  earthquakes  were 
found  in  Grand  Forks  AF'B.  The  difference  in  the  number  of  reported  events  is 
due  to  the  preponderance  of  earthquakes  in  the  600  to  1000  km  range  from  Minot, 
including  the  activity  from  the  Montana,  Yellowstone,  and  portions  of  the  Wvoming 
and  Colorado  source  regions  depicted  in  Figure  12.  The  seismic  sources  shown 
in  this  figure  were  used  to  evaluate  the  seismic  risk  at  Minot  AFB  and  portions  of 
the  last  two  were  within  the  1000  km  radius  used  to  evaluate  the  risk  at  Grand  Forks 
AFB.  The  seismicity  parameters  for  these  source  regions  are  given  in  Table  7. 

Outside  of  these  sources,  the  seismic  activity  within  1000  km  of  both  basics 
has  a  low  incidence  of  occurrence  with  a  wide  scatter  of  earthquake  epicenters. 
Together  these  eondi  lions  make  the  identification  of  seismic  trends  or  souree 
areas  impossible.  In  this  situation  il  is  neeessarv  to  assume  that  the  activitv 
occurs  randomly  throughout  the  entire  region  mtside  of  explicit  sources.  All 
activity  not  included  in  the  source  regions  listed  above  were  lumped  into  one  data 
set  for  evaluation  of  a  single  recurrence  <  ut  ve.  This  eurve  was  used  as  the 
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background  activity  level  for  all  sites  except  Whiteman  A  I  'M  .m  l  the  parameters 
of  the  curve  are  given  in  Table  7  under  the  label  Background. 

I  bis  assumption  of  random  distribution  of  activity  could  be  conservative  m 
the  sense  that  the  activity  is  assumed  to  have  an  equal  probal  :  1  i T v  of  occurring 
anvwhere  including  at  the  site  of  interest.  While  the  lack  of  sufl'i*  ient  knowledge 
concerning  the  earthquake  activity  in  this  region  precludes  any  other  approa>dt,  it 
is  possible  that  the  conditions  required  to  generate  an  earthquake  do  not  exist 
anywhere  near  the  sites  of  interest  and  would  produce  a  reduction  of  the  seismic 
hazard  at  these  wings. 

It  is  apparent  from  the  range  of  the  closest  seismic  source  region,  approxi¬ 
mately  500  km  for  Minot  .A  KB  and  750  km  for  Cirand  Forks  Al'U,  that  the  back¬ 
ground  seismicity  dominates  the  seismic  hazard  at  each  site  and  that  the  hazard 
is  essentially  equal  at  both  facilities.  The  difference  being  so  small  that  they  must 
he  considered  negligible.  In  Figure  17.  the  site  intensity  risk  curve  for  these 
liases  is  plotted  and  in  Figure  18  the  peak  ground  acceleration,  velocity,  and  dis¬ 
placement  curves  are  shown.  The  ground  motion  amplitudes  at  specific  levels  of 
annual  risk  are  also  given  in  I  able  10  and  in  Table  11  the  ground  motion  levels 
estimated  on  the  basis  of  intensity  are  given.  For  acceleration  the  values  arc 
within  a  factor  of  two  throughout  the  entire  range  and  can  be  considered  in  good 
agreement.  Both  velocity  and  displacement  tend  to  converge  at  the  lower  levels 
of  risk  but  are  substantially  different  at  the  higher  risks. 

In  Figures  10a,  19b,  anil  19c,  the  composite  horizontal  response  spectra  for 
10-,  100-,  and  1000-vear  return  period  ground  motions  are  shown.  These  figures 
are  based  on  the  ground  motions  given  in  Table  10. 

The  closest  reported  earthquake  activity  to  either  facility  was  an  epicomral 
intensity  IV  earthquake  reported  at  VVilliston,  .North  Dakota  in  26  October  1946  and 
was  located  approximately  17  5  km  to  the  west  of  Minot  A  FB,  nearly  on  the  Montana 
North  Dakota  state  line.  The  nearest  activity  to  Cirand  Forks  A  FB  occurred  on 
23  December  1028  in  north  central  Minnesota  at  an  epicentral  distance  of  260  km. 
Kach  site  is  within  300  km  of  several  other  events  all  of  which  W'cre  small  to 
moderate  in  magnitude  with  the  largest  event  having  an  epicentral  intensity  of  VII. 
As  stated  earlier,  the  locations  of  the  activity  is  erratic  and  provides  little  informa 
tion  tin  iho  location  of  future  activity  near  (he  facilities.  The  epicenters  of  the 
earthquakes  within  300  km  of  Minot  AI-'B  and  Cirand  Forks  AFB  arc  plotted  in 
Figures  20  and  21,  respectively.  These  events  cover  the  period  of  1917  to  1970. 
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Figure  18.  Peak  G round  Acceleration,  Velocity,  and  Displacement 
Seismic  Risk  Curves  for  Minot  and  Grand  Forks  Air  Force  Bases 


fable  10.  Peak  Ground  Motion  Annual  Risk  Levels  for  Mi  ,ot  A F H 
and  Grand  Forks  AFB 


- - 

Annual 

Risk 

Return 
Period 
(years ) 

Intensity 

Acceleration 
(cm  ' see- ) 

Velocity 
(cm  /  sec) 

Displacement 

(cm) 

0.  0 

1.  1  1 

— 

1.  5 

0.  06 

0.  02 

0.  5 

2 

I 

2.  4 

0.  14 

0.  04 

0.  2 

5 

11 

6.  4 

0.  28 

0.  08 

0.  1 

10 

II 

0.  3 

0.  42 

0.  13 

0.  05 

20 

III 

13.  1 

0.  61 

0.  18 

0.  02 

50 

III 

20.  3 

1.  00 

0.  29 

0.  0  1 

100 

III-IV 

27.  9 

1.  43 

0.  40 

0.  005 

200 

IV 

38.  2 

2.  06 

0.  54 

0.  002 

500 

I\  -\ 

57.  3 

3.  3  0 

0.  82 

0.  00  1 

1000 

\ 

77.  2 

4.  66 

1.11 
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Table  11.  Peak  Ground  Motion  .Annual  Risk  Levels  at  Minot  AFU 
and  Grand  Forks  AFU  Based  on  Intensity  Levels 


Annual 

Risk 

Fractional 

Intensity 

Acceleration 
(cm  /  sec^ ) 

Velocity 
(cm  /  sec) 

Displacement 

(cm) 

0.  9 

— 

— 

1  1 

--- 

0.  5 

1.  4 

2.  63 

0.  17 

0.  2 

2.  2 

4.  55 

wtm 

0.  26 

0.  1 

2.  7 

6.  18 

1.  10 

0.  33 

0.  05 

3.  1 

8.  27 

1.  40 

0.  4  1 

0.  02 

3.  6 

11.  56 

1.  85 

0.  54 

0.  01 

4.  0 

13.  93 

2.  16 

0.  63 

0.  005 

4.  3 

19.  14 

2.  82 

0.  81 

0.  002 

4.  8 

26.  85 

3.74 

1. 06 

0.  001 

5.  2 

34.20 

4.  57 

1.  28 

Figure  19a.  Composite  Horizontal  Design  Response 
Spectra  for  Minot  and  Grand  Forks  Air  Force  Bases 
Using  10 -year  Return  Period  Ground  Motions 
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Figure  19b.  Composite  Horizontal  Design  Response 
Spectra  for  Minot  and  Grand  Forks  Air  Force  Bases 
Using  100-year  Return  Period  Ground  Motions 


Figure  20.  Reported  Earthquake  Epicenters 
Within  300  km  of  Minot  AFB 
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Figure  21.  Reported  Earthquake  Epicenters 
Within  300  km  of  Grand  Forks  AFB 


3.4  Whiteman  AFB  -  Wing  IV 


Whiteman  AFB  is  located  near  what  is  probably  the  most  active  seismic  region 
in  the  United  States  east  of  the  Rocky  Mountains.  The  area  near  the  confluences 
of  the  Ohio.  Missouri,  and  Mississippi  Rivers  encompasses  a  large  percentage  of 
the  epicenters  of  earthquakes  reported  in  the  central  United  Slates.  Among  the 
events  reported  in  this  region  are  several  of  the  largest  events  to  have  occurred 
in  the  continental  United  States  during  historic  times.  The  location  of  Whiteman 
AFB  relative  to  this  source  region,  the  Central  Mississippi  Valley  area,  is  shown 
in  Figure  22.  In  addition  to  this  source  region,  only  a  source  region  designated 
as  the  Southeastern  United  States  has  any  extent  within  the  1000  km  radius  of 
Whiteman  AFB.  Outside  of  these  two  regions  the  seismic  activity  in  the  central 
United  States  appears  to  be  randomly  distributed  and  was  treated  as  background 
seismicity.  All  risk  calculations  were  conducted  for  38.  96°N  and  92.  35°W,  the 
coordinates  of  Columbia,  Missouri.  As  discussed  at  the  end  of  Section  2.  3,  some 
modification  of  the  risk  estimate  for  this  site  is  expected  due  to  changes  presently 
being  made  to  the  attenuation  function  modification  procedure. 


Figure  22.  Seismic  Source  Regions  Used  in 
the  Evaluation  of  Seismic  Risk  at  Whiteman 
AFB  (indicated  by  closed  circle) 


Previous  seismicity  studies  have  been  conducted  for  both  of  the  source  regions 
included  in  risk  evaluations  for  Whiteman  AFB  and  these  studies  were  used  to 
define  the  recurrence  curves  for  the  source  regions.  The  Southeastern  United 
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States  source  area  was  investigated  by  Bollinger.  The  recurrence  curve 

generated  in  this  study  was  converted  by  Modified  Mercalli  intensity  to  local  mag- 

14 

nitude,  M  using  a  relation  derived  for  the  eastern  United  States  by  Brazee. 

For  the  Central  Mississippi  Valley  region  several  studies  have  been  conducted  and 

i,  15.16.17 

the  recurrence  curves  are  in  reasonable  agreement.  1  he  recurrence 

17 

curves  developed  by  Battis  for  this  source  region  and  for  the  central  United 
States  background  seismieitv  were  used  in  this  report.  The  parameters  defining 
the  seismic  activity  for  these  areas  are  given  in  Table  12. 


Table  12.  Whiteman  APB  Source  Region  Parameters 


Source 

Area 
<104  km) 

log  (N/V)  =  A  -  bAJ  ^ 

,,  MAX 

A  b 

Central  Mississippi  Valley 

15.  9 

2.  12  0.  62 

8.  0 

Southeastern  United  States 

43.  2 

3.09  0.85 

7.  5 

Background 

— 

-0.  58  ^  0.56 

6.  0 

^ log  ( N / y  / 1 04  km2) 


The  seismic  hazard  for  Whiteman  AFB  was  calculated  using  both  site  intensity 
and  regionally  modified  acceleration  attenuation  functions  for  the  central  United 
States  as  given  in  Table  2.  At  present,  no  satisfactory  independent  means  is 
available  to  estimate  peak  ground  velocity  or  oispiacement  risk  curves  for  the 
central  United  States  and  these  curves  could  only  be  obtained  by  scaling  of  either 
the  acceleration  or  intensity  risk  curves.  Scaling  based  on  intensity  is  done  using 
Kqs.  (3),  (4),  and  (5).  Given  the  acceleration  risk  curve,  velocity  and  displacement 
curves  can  be  estimated  using  the  "standard"  earthquake  values  as  determined  by 
Newmark  and  Hall.  The  "standard"  earthquake  assumes  a  ratio  of  acceleration 
to  velocity  and  displacement  of  0.5:60.96  cm/sec:45.72  cm.  The  same  t  atios  are 
assumed  for  all  other  levels  of  acceleration.  This  procedure  is  open  to  question 
for  several  reasons.  Those  ratios  were  determined  by  analysis  of  records  of  the 
1940  Imperial  Valley,  California  earthquake.  The  ratio  of  high  frequency,  motions, 
or  acceleration,  to  longer  period  motions  of  velocity  and  displacement  at  the  source 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  57.  ) 
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appear  to  be  dependent  on  many  factors  such  as  fault  rupture  length,  rise  time, 
and  stress  drop  which  vary  with  the  magnitude  of  the  earthquake  and  regional 
tectonics.  Away  from  the  source,  regional  variation  in  frequency  dependent 
attenuation  will  also  affect  the  validity  of  these  ratios.  However,  there  does  not 
appear  to  be  a  clearly  superior  method  for  evaluating  the  velocity  and  displace¬ 
ment  risk  curves. 

The  site  intensity  risk  curves  and  evaluated  peak  acceleration  and  scaled 
velocity  and  displacement  risk  curves  are  shown,  as  the  curves  labeled  "a",  in 
Figures  23  and  24  respectively.  Values  of  these  ground  motion  descriptors  at 
specified  levels  of  annual  risk  are  also  given  in  Table  13.  In  Table  14,  the  equiva¬ 
lent  peak  ground  motion  values  determined  by  conversion  from  site  intensity  are 
given.  Comparison  of  the  two  sets  of  values  given  in  Table  14  indicate  a  large 
discrepancy  between  the  seismic  hazard  as  indicated  by  site  intensity  and  peak 
ground  acceleration.  To  determine  which  risk  curve  is  more  reasonable,  a  com¬ 
parison  was  made  between  the  predicted  recurrence  periods  for  site  intensities  and 

the  periods  based  on  historic  observations.  The  source  for  the  historic  data  was 

14 

a  compilation  by  Brazee  covering  a  one -degree  square,  including  Columbia,  for 
the  45-year  period  from  1928  to  1973. 

Historically,  the  Columbia  area  has  experienced  intensity  Ill  or  greater  six 
times,  IV  or  greater  on  four  occasions,  two  V  or  greater,  and  one  intensity  \'I 
during  this  period.  The  intensity  based  risk  curve  predicts  24  occurrences  of  111 
or  greater,  12  of  IV  or  larger,  6  V  or  greater,  2  intensity  VI  or  greater  and 
1  intensity  VII  or  larger  earthquake.  Though  the  actual  historic  anti  predicted 
rates  should  not  be  expected  to  parallel  exactly,  the  departure  between  historic 
record  and  prediction  is  significant  especially  at  the  low  intensities  which,  due  to 
the  frequency  of  occurrence  should  be  reflected  more  accurately  in  the  risk  curve. 
Using  Kq.  (3)  the  mean  value  of  acceleration  at  each  intensity  can  be  determined 
and  the  expected  number  of  occurrences  in  the  45-year  period  determined.  These 
calculations  yield  9,  4,  2  and  I  occurrences  of  intensities  III,  IV,  V  and  VI  or 
greater,  respectively,  at  Whiteman  AFB.  While  this  comparison  is  only  approxi¬ 
mate,  it  suggests  that  the  intensity  risk  curve  overstates  the  risk  at  Whiteman  AFB 
and  the  acceleration  risk  curve  appears  to  be  more  realistic  in  representing  the 
ground  motions  experienced  at  this  facility. 
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MODIFIED  MERCALLI  INTENSITY 


Figure  23.  Modified  Mercatli  Intensity  Seismic  disk  Oirv< 
lor  Whiteman  A  KB.  Curve  (a)  assumes  the  maximum 
magnitude  earthquake  can  occur  anvwherc  in  the  Central 
Mississippi  Valiev  source  region.  Curve  (h)  restricts  this 
event  to  the  southern  extent  of  the  source  region.  See  text 
for  fuller  discussion. 


Table  14.  Peak  Ground  Motion  Annual  Risk  Levels  for 
Whiteman  AFB  Based  on  Intensities 


r.  _ _ _ 

Annual 

Risk 

Fractional 

Intensity 

Acceleration 

(cm/sec^) 

Velocity 
(cm /sett) 

I  iisplacement 
(cm ) 

0.  9 

0.  5 

2.  6 

5.  8 

1.  0 

0.  3 

0.  2 

4.  2 

17.7 

2.  6 

0.  8 

0.  1 

5.  2 

35.  2 

4.7 

1.  3 

0.  05 

6.  1 

65.  0 

7.  8 

2.  1 

0.  02 

T.  3 

152.  1 

15.  8 

4.  2 

0.  01 

8.  2 

27  9.  3 

26.  3 

6.  ! * 

0.  005 

!>.  0 

498.  9 

42.  7 

11.0 

0.  002 

10.  0 

968.  0 

74.  1 

18.  6 

0.  001 

10.  53 

1396.  4 

100.  6 

25.  0 

As  with  most  of  the  areas  studied  in  this  report  the  direct  correlation  of 
earthquake  epicenters  with  known  faults  in  the  Central  Mississippi  Valiev  region 
is  not  possible  at  the  present  time.  In  this  sense,  the  boundaries  of  source  region 
are  arbitrary  in  their  assignment.  To  indicate  the  possible  variation  in  the 
seismic  hazard  at  Whiteman  AFB  due  to  this  factor  a  contour  plot  of  the  500 -year 
return  period  accelerations  for  a  one -degree  square  centered  on  Whiteman  AFB 
is  shown  in  Figure  25.  By  reeiprosity,  expansion  of  the'  source  region  towards  the 
base  is  equivalent  to  moving  the  base  towards  the  source.  Conversely,  reduction 
of  the  source  area  is  equivalent  to  moving  the  site  of  interest  to  the'  northwest  of 
its  present  location.  The  gradient  in  this  region  is  not  very  steep  and  changes  of 
the  order  of  a  few  tens  of  kilometers  in  the  bounds  of  the  Central  Mississippi 
Valley  cource  region  would  not  substantially  change  the  seismic  hazard  estimate 
at  Whiteman  AFB. 

A  second  potential  modifier  of  the  risk  curves  involves  the  non-uniformity  of 
earthquake  activity  within  the  Central  Mississppi  Valley  source  region.  Insufficient 
data  is  presently  available  to  subdivide  this  region  into  smaller  areas,  and  thus 
account  for  the  non -uniformity  of  the  activity.  However,  another  approach  was 
used  to  obtain  a  limiting  lower  bound  to  the  risk  curves.  It  was  assumed  that  the 
earthquake  process  in  much  of  the  Central  Mississippi  Valiev  source  region  could 
not  generate  the  maximum  magnitude  earthquake  of  !i.  0  Mj  for  this  region.  These 
large  events  were  restricted  to  a  region  at  the  southern  end  of  the  source  area  in 
the  vicinity  of  the  Missouri  -  Arkansas  -  Tennessee  junction.  For  the  remainder 
of  the  source  area  a  maximum  magnitude  earthquake  of  6.5  Mj  was  allowed.  The 
seismic  risk  curves  for  this  case  are  given  in  Figures  23  and  24  as  the  curves 


labeled  "b".  In  Table  15,  the  ground  motion  levels  at  specified  annual  risks  are 
given.  These  should  be  considered  only  as  the  lower  limiting  risk  with  the  true 
curve  lying  somewhere  between  these  two  curves. 


92.5*  92*  W 


Figure  25.  Contours  of  the  500-year 
Return  Period  Accelerations  in  a 
One -degree  Square  centered  on 
Whiteman  AFH  (closed  circle); 

(unit.;  of  cm/sec^).  Accelerations 
calculated  assuming  the  maximum 
magnitude  earthquake  can  occur 
anywhere  within  the  Central 
Mississippi  Valley  source  region 


Table  15.  Lower  Limit  Peak  Ground  Motion  Annual  Risk  Levels  of 
Whiteman  A FB 


i  -  — - - : - - - 1 

Annual 

Return 

Period 

Intensity 

Acceleration 

(om/siv2) 

V  elocitv 
(cm  /sec) 

Displacement 
(cm ) 

0.  9 

1.  11 

0.  8 

0.  1 

0.  1 

0.  5 

2 

11 

2.  5 

0.  3 

0.  2 

0.  2 

5 

IV 

6.  9 

0.  9 

0.7 

0.  1 

10 

V 

13.  0 

I.  6 

1.  2 

0.05 

20 

V  -VI 

22.  7 

2.  8 

2.  1 

0.  02 

50 

VI -VII 

43.  8 

5.  5 

4.  1 

0.  01 

100 

V  11 

67. 7 

8.  4 

6.  3 

0.  005 

200 

VII  -VIII 

100.  7 

12.  6 

4 

0.  002 

500 

VIII 

161.5 

20.  1 

15.  1 

0.  001 

1000 

VIII-IX 

224. 7 

28.  0 

2  1.0 

The  composite  horizontal  design  response  spectra  lor  1(1  ,  1(H)-,  and 
1000 -vear  return  periods  are  shown  in  Figures  26a,  261»,  and  26c.  These  spectra 
are  based  on  the  ground  motions  given  in  Table  lit  which  assumed  a  uniform  potential 
for  a  magnitude  8.  0  M  .  earthquake  throughout  Cent  ml  M  ississippi  \  a  I  lev  Source  Region 
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Figure*  26a.  Horizontal 
Design  Response  Spectra 
for  Whiteman  AFB  Based 
on  10-vear  Return  Period 
Ground  Motions  Assuming 
the  Maximum  Magnitude 
Harthquake  Can  Occur 
Anywhere  in  the  Central 
Mississippi  Valiev  Source 
Region 
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Figure  26b.  Horizontal 
Design  Response  Spectr  i 
for  Whiteman  AFB  Based 
>n  100 -year  Return  period 
Ground  Motions  Assuming 
the  Maximum  Magnitude 
Karthquake  Can  Occur 
Anywhere  in  the  Central 
Mississippi  \  allev  Sour-  e 
Region 
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Figure  26r.  Horizontal 
Design  Response  Spectra 
for  Whiteman  AFB  Based 
on  10 00 -year  Return  Period 
Ground  Motions  Assuming 
the  Maximum  Magnitude 
earthquake  Gan  Occur 
Anvwhere  in  tlie  Central 
Mississippi  Valiev  Source 
Region 


In  Figure  27,  all  earthquakes  reported  in  the  Karthquake  Data  File  with  epi¬ 
centers  within  200  km  of  Whiteman  AFB  are  plotted.  Of  these  27  events,  the 
largest  reported  magnitude  is  one  of  5.  6  which  occurred  in  March  1974  at  a 
distance  of  198  km  from  the  base.  The  nearest  event  is  at  108  km  but  has  no 
reported  magnitude  or  intensity.  The  concentration  of  epicenters  to  the  south¬ 
east  of  Whiteman  AFB  lays  within  the  Central  Mississippi  Valley  seismic  source 
/•one.  Several  earthquakes  not  shown  in  this  figure  should  also  be  noted.  The 
most  important  of  these  was  the  1811-1812  series  of  shocks  which  occurred  near 
New  Madrid,  Missouri  in  southeastern  Missouri  and  northeastern  Arkansas. 

Nuttli  estimated  the  magnitude  of  the  three  principal  shocks  to  be  7.  3,  7.  1  and 

19 

7.4  m.  .  1  In  the  meizoseismal  zone  it  is  thought  that  Modified  Mercalli  intensity 

XII  was  reached.  The  first  of  these  events  was  reported  felt  as  far  away  as  Boston 

20 

at  an  epicentral  distance  of  approximately  1800  km.  .  For  purposes  of  compari¬ 
son,  the  area  enclosed  in  the  Modified  Mercalli  intensity  V  isoseismal  was  esti- 

6  2 

mated  by  Nuttli  to  be  2,  5  X  10  km  while  the  1960  San  Francisco  earthquake  pro- 

5  2  19 

duced  this  intensity  over  an  area  of  only  1.5  X  10  km  .  Assuming  that  the 

maximum  magnitude  earthquake,  8.0  M(  ,  were  to  reoccur  near  the  epicenter  of 

these  events,  the  predicted  ground  motions  at  Whiteman  AFB  would  be  approxi - 

mately  428  cm/sec",  53  cm/sec  and  40  cm.  This  corresponds  to  a  site  intensity 

of  VIII  -  IX,  approximately  one  unit  higher  than  that  obtained  by  extrapolation  of 

1 9 

Nuttli 1  s  isoseismal  map  for  the  first  event  in  this  series.  "  The  horizontal  design 
response  spectra  for  this  event  does  not  exceed  the  1000-year  return  period 
spectra  shown  in  Figure  26c. 

3.5  K.K.  Warren  AFB  -  Wing  V 

The  seismic  source  regions  used  in  the  hazards  evaluation  for  F.  K.  Warren 
AFB  are  identical  to  those'  used  for  FUsworth  AFB  and  are  shown  in  Figure  12. 
These'  regions  are  based  on  an  examination  of  all  seismic-  activity  within  1000  km 
of  the  facility.  The  associated  seismicity  parameters  for  these  areas  are  given  in 
Table  7.  The  coordinates  of  Cheyenne,  Wyoming,  4  1.  1°N  and  104.  9° W  were  used 
as  the  location  for  which  the  seismic  risk  was  calculated. 

The  annual  risk  curves  for  F.  K.  Warren  AFB,  for  Modified  Mercalli  intensity 
and  peak  ground  acceleration,  velocity,  and  displacement  are  shown  in  Figures  28 
and  29,  respectively.  As  with  all  sites  except  for  Whiteman  AFB,  each  risk  curve 
is  independently  evaluated  using  the  appropriate  attenuation  function  as  specified 


lit.  Nuttli,  ().  (  1973)  The  Mississippi  Valley  earthquakes  of  1811  and  1812; 
Intensities,  ground  motion  and  magnitudes.  Bull.  Seism.  Soc.  Am. 

63:227  -24H.  “  ~  ' 

20.  V  nn  Hake,  (',  (1974)  Karthquake  history  of  Missouri ,  Karthquake  Information 
Bull.  6:25-26.  ' 
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Figure  27.  Reported  Farthquake  Kpieenters 
Within  200  km  of  Whiteman  AFB 

in  Table  2.  For  specific  levels  of  annual  risk,  the  calculated  ground  motions  at 
F.  F.  Warren  Al-'B  are  given  in  Table  13.  In  Table  17  the  ground  motions  for  these 
same  annual  risks,  as  evaluated  by  conversion  of  site  intensities,  are  given.  The 
two  sets  of  values,  given  the  limitations  of  the  conversion  procedure,  are  considered 
to  be  in  reasonable  agreement.  The  ground  motions  calculated  directly,  however, 
are  the  preferred  values. 

As  can  be  seen  in  Figure  12  the  location  of  K.  K.  Warren  Al-'B  is  in  close  proxi  - 
metrv  to  the  boundaries  of  two  source  regions,  those  designated  as  Colorado  and 
Wyoming  in  this  figure.  In  Figure  30  the  contours  of  the  300-vear  year  return 
period  accelerations  within  a  one -degree  square  centered  on  Wing  Y  are  show  n. 

The  general  characteristics  of  this  field  also  hold  for  other  risk  levels  and  ground 
motion  parameters  though  typically  the  gradient  is  shallower.  Beyond  the  spatial 
variation  of  seismic  risk  in  the  vicinity  of  F.  F.  Warren  AFB,  this  figure  also 
provides  information  on  the  effects  of  relocation  of  the  borders  of  the  Colorado 
and  Wyoming  source  regions.  Movement  of  the  site  towards  or  away  from  tin' 
boundary  of  these  source  regions  is  equivalent  to  moving  the  borders  relative  to 
the  site.  As  these  borders  are  arbitrarily  assigned  due  to  various  factors  dis¬ 
cussed  in  Section  2  of  this  report,  this  figure  provides  an  indication  of  the  uncer¬ 
tainty  in  the  risk  evaluation  which  is  associated  with  these  factors. 
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Kitfure  28.  Modified  Morcalli  Intensity  Seismic  Risk  Curve  for 
K.  F.  Warren  A  KM 


PEAK  GROUND  MOTION 


Figure  2!'.  Peak  (i round  Acceleration,  Velocity,  and  Displacement  Seismic 
Risk  Curves  for  F.  K.  Warren  AFB 


Table  16.  Peak  Ground  Motion  Annual  Hi.sk  Levels  for  F.  K.  Warren  AFB 


Annual 

Risk 

Return 
Period 
(years ) 

Intensity 

Acceleration 
(cm  /  soc^ ) 

Velocity 
(cm  /  sec ) 

Displacement 

(cm) 

0.  0 

1.  11 

l-li 

6.  7 

0.  2 

0.  05 

0.  5 

> 

11 

10.  8 

0.  5 

0.  12 

0.  2 

- 

HI 

1G.  1 

o.  n 

0.  26 

0.  1 

10 

!\ 

21.0 

1.  4 

0.  -1  I 

0.  05 

20 

IV 

27.  2 

2.  0 

0.  61 

0.  02 

50 

IV  -V 

28.  5 

2.  2 

0.  07 

0.  0  I 

100 

V 

50.  1 

■1.  5 

1.  25 

0.  005 

200 

\ 

65.  2 

6.  2 

1.  84 

0,  002 

500 

V  1 

02.  1 

!».  :■{ 

2.70 

0.  001 

1000 

V  I 

121.8 

12.  4 

2.  55 

Table  17.  Peak  Ground  Motion  Annual  Risk  Levels  for 
F.  E.  Warren  AFB  Based  on  Intensity  Levels 


Annual 

Risk 

F  factional 
Intensitv 

Acceleration 

/set.'2) 

Velocity 
(cm  /  sec) 

Displacement 
(cm ) 

0.  9 

1.  8 

.1.  3 

0.7 

0.  2 

0.  5 

2.  7 

6.  :i 

1.  1 

0.  3 

0.  2 

3.  5 

li.  l 

1.  8 

0.  5 

0.  1 

4.  0 

15.  7 

2.  4 

0.  7 

0.  05 

4.  4 

20.  7 

3.  0 

0.  9 

0.  02 

5.  0 

29.  9 

4.  1 

1.  2 

0.  01 

5.  3 

37.7 

5.  0 

1.  4 

0.  005 

5.  6 

47.  4 

6.0 

1.  7 

0.  002 

6.  1 

7.7 

2.  1 

0.  001 

6.  4 

9.  1 

2.  5 
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Figure  BO.  Contours  of  the 
500-year  Return  Period 
Acceleration  in  a  One -degree 
Square  Centered  on 
F.  K.  Warren  AFB  (closed  circle). 
41.0'  (f’nits  of  cm/sec2) 


In  Figures  .'11a,  .'lib,  and  tile  the  composite  horizontal  design  response  site,  fa 
for  10-,  100-,  and  1000-year  return  period  ground  motions  as  given  in  '1  able  la, 
are  sho.t  . 


VELOCITY  CCM/5EC)  VELOCITY  ( CM/5EC  )  VELOCITY  CCM/5EC) 


Figure  3la.  Composite 
Horizontal  Design 
Response  Spectra  lor 
F.  E.  Warren  AFH  Based 
on  10 -year  Retur  n 
Period  Ground  Motions 
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Figure  31b.  Composite 
Horizontal  Design 
Response  Spectra  lor 
F.  E.  Warren  AFB  Based 
on  100 -year  Return 
Period  Ground  Motions 


Figure  31c.  Composite 
Horizontal  Design 
Response  Spectra  lor 
F.  E.  Warren  AFB  Based 
on  1000 -year  Return 
Period  Ground  Motions 


EREHLENCY  (H?) 


Figure  32  is  a  plot  of  all  earthquakes  which  have  occurred  within  200  km  of 
F.  E.  Warren  AFH.  A  total  of  63  events  arc*  plotted  covering  a  period  from  1692 
to  107  8.  The  closest  event  was  located  approximately  8  km  from  the  site  of 
interest  but  had  no  reported  magnitude  or  intensity.  The  largest  event  is  one  of 
magnitude  5.3  m,|  and  epicentral  intensity  of  \  II  which  .m  urr.-d  jn  August  10C7. 
It  was  located  at  135  km  front  K.  E.  Warren  A  KM  in  tim  cluster  of  epicenters  due 
south  of  the  site.  The  vast  majority  of  the  events  reported  have  body-wave  mag¬ 
nitudes  between  4.  0  and  5.3  and  were  reported  since  1963.  Onlv  five  earth¬ 
quakes  are  reported  prior  to  1963  indicating  the  incompleteness  of  the  data  set. 
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Figure  32.  Reported  Earthquake  Epicenters 
Within  200  km  of  E.  E.  Warren  AFH 


\  irtually  all  of  the  activity  located  in  the  cluster  due  south  of  the  site  is 
associated  with  fluid  injection  at  the  Rocky  Mountain  Arsenal  outside  of  Denver, 
Colorado.  It  is  assumed  that  the  deep  well  injection  of  waste  fluids  from  the 
Arsenal  operation  acted  as  a  lubricant  on  faults  within  this  region  and  allowed  the 
rupturing  of  faults  which  would  otherwise  have  remained  locked  until  greater 
stresses  could  have  accumulated.  In  a  sense,  this  indicates  that  the  inclusion  of 
area  in  the  Colorado  source  region  may  be  inappropriate  as  these  appear  to  be 
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induced  earthquakes.  However,  at  least  two  earthquakes  have  occurred  in  this 
region  prior  to  the  fluid  injections,  in  1882  and  1924,  and  the  post-injection 
activity  itself  suggests  that  accumulation  of  sufficient  strain  in  this  region  for  sig¬ 
nificant  future  seismic  activity.  Thus,  it  is  felt  that  this  area  should  be  included 
in  the  Colorado  source  region. 

4.  CONCLUSIONS 

Using  standard  methods  of  probabi listic  seismic  risk  analysis,  estimates  of 
the  seismic  hazard  at  six  Minuteman  missile  wings  were  conducted.  For  four  of 
the  six  sites,  the  dominant  contribution  to  the  site  hazard  is  the  general  background 
•seismicity.  These  sites  were  Malstrom,  KUsworth,  Minot,  and  Grand  Forks  Air 
Force  liases.  The  background  activity  is  used  to  represent  earthquakes  which 
appear  to  occur  randomly  throughout  the  regions  surrounding  the  sites  and  which 
can  not  be  associated  with  any  distinct  source  region.  The  apparent  randomness 
of  these  events  is,  at  Least  in  part,  the  result  of  low  rates  of  occurrence  and  a 
short  historical  record  which,  taken  together,  greatly  limit  the  ability  of  the  in¬ 
vestigator  to  isolate  these  events  within  a  source  region.  The  assumption  of 
random  occurrence  for  these  events  which  was  used  in  these  studies,  is  to  some 
degree,  conservative  aS  this  activity  may  be  more  limited  in  areal  extent  and  thus 
the  estimated  risk  curves  may  also  be  conservative. 

The  site  having  the  greatest  seismic  hazard  is  Whiteman  AFB,  Missouri. 
Though  not  commonly  considered  to  be  a  region  of  high  seismic  hazard  due  to  the 
lack  of  recent  major  earthquakes,  the  combination  of  low  seismic  wave  attenuation 
in  the  central  United  States  and  the  potential  for  major  earthquakes  in  the  central 
Mississippi  Valley  combine  to  produce  high  expected  ground  motions  at  this  site  for 
long  inter -occurrence  periods.  For  short  return  periods,  less  than  20  years,  the 
expected  ground  motions  are  lower  than  those  at  the  western  United  States  missile 
wings.  The  uncertainties  in  ground  motion  estimates  for  this  region  are  larger 
than  for  the  vest  coast  and  it  is  expected  that  future  studies  will  modify  the  results 
presented  in  this  report.  The  changes  presently  being  incorporated  into  the  tech¬ 
nique  used  to  generate  the  peak  acceleration  attenuation  function  which  was  used 
for  this  site  are  expected  to  reduce  the  seismic  hazard  although  a  quantitative 
evaluation  of  the  effect  is  not  possible  at  present. 

F.F.  Warren  AFB  has  the  highest  level  of  risk  at  the  remaining  sites.  Both 
for  this  wing  and  for  Malstrom  AFB,  significant  changes  in  the  predicted  levels  of 
risk  can  be  obtained  bv  relatively  minor  changes  in  the  boundaries  of  close  seismic 
source  regions.  As  t he*  boundaries  of  those  regions  are  not  well  established,  the 
variability  associated  with  alteration  of  the  boundaries  must  be  ronsulrm)  in  an v 
application  of  these  risk  estimates. 
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Appendix  A 

Regional  Modification  of  Acceleration 
Attenuation  Functions 

Al.  INTRODUCTION 

To  rrmvcrt  the  spatial  and  temporal  charaeterist  i  i  s  .>)'  seismic  activity  into  an 
estimate  of  the  seismic  hazard  at  the  site  of  interest  it  is  neecssarv  to  have  a 
method  for  the  determination  of  ground  motion  at  a  location  remote  from  the  epi¬ 
center  of  an  earthquake  of  specified  magnitude.  A  gem  ral  set  of  functions  of  : In¬ 
form: 

g  GOU.H)  (All 

■^s 

where  g^  is  the  ground  motion  al  the  site  of  interest,  M  is  tin-  earthquake  magni¬ 
tude  and  K  is  the  site  epieentral  or  hvpcic-entr.il  distance,  have  heen  derived 
empiric-allv  for  this  purpose.  Typically,  the  ground  motion  descriptor,  g ,  is  m 

terms  of  peak  acceleration,  velocity,  or  displacement;  though  other  measures 
have  been  used  (McGuire,  107 6 f . 

A  problem  if  great  significance  is  associated  with  the  use  of  these  functions. 
Most  of  the  data  upon  which  these  equations  tire  based  are  obtained  in  California. 
Onlv  a  very  limited  number  of  strong  motion  records  tire  available  from  the  rest 
of  the  continental  United  States.  At  the  same  time,  there  ire  manv  lines  of  evi¬ 
dence  that  ground  motion  attenuation  in  California  is  atypical  compared  to  other 
Note:  1-iee  page  <2  for  Appendix  A  References. 
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areas  of  the  United  States  (Brazee,  1976;  NuttLi,  1 97 3 ) .  Thus,  the  direct  applica¬ 
tion  of  attenuation  functions  based  primarily  on  western  United  States  observations 
to  other  regions  is  at  best  only  a  gross  approximation  of  the  real  situation. 


A2.  PAST  MODIFICATION  SCHEMF.S 

Several  methods  have  been  suggested  to  overcome  the  difficulties  of  regional 
variation  in  ground  motion  attenuation.  The  first  solution  was  the  use  of  Modified 
Mercalli  intensity  as  the  ground  motion  descriptor,  g  .  Intensity  attenuation 
functions  of  the  form; 

I  =  I(M,  r)  or  1(1  ,  K)  (A2) 

S  o 

where  I  is  the  site  intensity  and  I()  the  epicentral  intensity,  replaced  the  ground 
niotion  attenuation  functions.  The  advantage  of  this  method  is  that  historical 
reports  of  earthquake  damage  can  be  interpreted  in  terms  of  Modified  Mercalli 
intensity,  greatly  expanding  the  data  base  for  analysis  as  compared  to  the  use  of 
ground  niotion  descriptors  requiring  instrumental  observations.  This  feature  is 
significant  for  areas  of  low  seismicity  such  as  the  central  or  eastern  United  States. 
In  addition,  the  use  of  intensity  appears  to  eliminate  the  dependence  on  western 
United  States  data  to  estimate  seismic  hazard. 

This  method,  however,  is  highly  unsatisfactory.  The  primary  problem  with 
this  technique  is  that  Modified  Mercalli  intensity  describes  ground  motion  on  a 
verv  coarse  scale  which  is  based  on  subjective  evaluations  of  damage.  Taken  alone, 
intensity  is  extremely  difficult  to  interpret  in  terms  of  the  engineering  implications. 
What  knowledge  that  does  exist  of  the  engineering  parameters  associated  with 
given  levels  of  Modified  Mercalli  intensity  is,  to  a  large  degree,  based  on  studies 
of  western  United  States  earthquakes.  Thus,  the  method  does  not  eliminate  the 
cross -regional  comparisons. 

To  eliminate  the  engineering  application  problems  of  using  intensity,  attempts 
have  been  made  to  correlate  site  intensity  directly  with  peak  acceleration,  velocity, 
arid  displacement  (Guttcnbcrg  and  Richter,  1942;  Trifunac  and  Brady,  1975; 

McGuire,  1!I77  ).  Given  a  correlation  between  one  of  the  ground  motion  parameters 
and  site  intensity; 

g  l'(I  )  (Alt) 

s  s 

it  would  be  possible  to  convert  the  intensity  attenuation  functions  to  ground  motion 
attenuation  as; 


60 


(A4) 


gs  =  K»  . 

It  is  obvious  from  the  data  shown  in  Figure  A1  that  the  expected  trends,  in¬ 
creasing  ground  motion  amplitude  with  increasing  site  intensity,  occur  on  the 
average.  However,  the  wide  spread  of  values  observed  for  a  given  intensity  level 
result  in  only  weak  correlations  between  each  of  the  ground  motion  parameters 
and  intensity  (Trifunac  and  Brady,  1975;  McGuire,  1977;  and  Murphy  and  O'Brien, 
1977).  Velocity  does  correlate  with  intensity  to  a  higher  degree  than  acceleration 
or  displacement.  Unfortunately,  most  seismic  engineering  codes  are  written  in 
terms  of  acceleration 
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Figure  Al.  Means  of  Observed 
Peak  Ground  Acceleration, 
Velocity,  and  Displacements, 
With  One  Standard  Deviation 
Krror  Bars,  Plotted  Against 
Modified  Mercalli  Intensity 


In  part,  the  low  correlation  can  be  explained  bv  the  coarseness  of  I  ho  Modified 
Mercalli  intensity  scale.  Another  significant  element,  however,  is  the  dispersion 
of  the  wave  train  as  it  travels  away  from  the  epicenter.  The  level  of  damage  which 
results  from  a  single,  impulsive  motion  can  also  result  from  a  lower  amplitude 
motion  acting  over  more  cycles.  To  incorporate  this  distance  effect,  correlations 
have  also  been  made  using  both  site  intensity  and  epieentral  distance; 


H  I 


(A5) 


g  G(I  ,  H>  . 

BS  S 

The  use  of  distance  significantly  improves  the  correlation  with  peak  acceleration 
and,  to  a  lesser  degree,  with  peak  displacement  (McGuire,  1977 ;  Murphy  and 
O'Brien,  1977). 

Still,  the  use  of  either  type  of  site  intensity-ground  motion  relationship,  Kqs. 
(A3)  or  (Ad),  outside  of  California,  in  conjunction  with  a'regional  intensity  attenua¬ 
tion  function,  provides  only  limited  improvement  over  the  direct  application  of 
California  based  ground  motion  attenuation  functions.  In  both  cases,  the  effect  ol 
distance  has  not  been  completely  removed.  In  the  first  form,  the  distance  effects 
are  completely  ignored.  In  the  second  type  of  relationship,  an  assumption  is  made 
that  the  dispersive  effects  are  regionally  invariant. 


AH.  THU  MKTHOI) 


To  overcome  the  problems  which  effect  intensity  based  ground  motion  esti¬ 
mates,  a  new  'echnique  is  proposed  for  the  estimation  of  acceleration  attenuation 
functions  for  regions  with  inadequate  strong  ground  motion  data  for  the  generation 
of  empirical  functions.  This  method  incorporates  a  site  intensity -acceleration 
correlation  but  on  a  much  more  restricted  basis  than  previously  used  techniques. 
A  limited  amount  of  instrumental  observations  is  required  but  in  a  form  which  can 
generally  be  obtained  for  regions  of  low  seismic  activity  and  short  histories  of 
instrumental  coverage.  The  method  is  dependent  on  three  assumptions. 

The  first,  and  least  significant  of  the  assumptions,  is  that  some  functional 
form  exists  which  can  depict  the  characteristics  of  peak  acceleration  attenuation 
in  any  region.  The  most  commonly  used  function  for  acceleration  attenuation 
equations  is  of  the  form: 


In  a  (V.  •  <v,.M  -  <».,  In  (It  *  It  ) 
s  1  2  i  o 


( A  6 ) 


where  a  is  the  site  peak  acceleration,  M  is  the  magnitude  of  the  causilative  earth¬ 
quake,  R  is  the  epicentral  distance,  R^  is  a  distance  factor  tvpicallv  set  to  either 

zero  or  25  km  and  <v  ,  <t ,  O'  are  regression  parameters  based  on  observational 
l  2  •> 

data.  In  this  report,  only  this  form  of  attenuation  lunction  will  be  considered,  1  he 
explicit  form,  however,  is  not  critical  to  the  method.  It  is  assumed  that  this 
equation,  given  the  appropriate  regression  parameters,  will  adequately  character¬ 
ize  attenuation  in  anv  region.  To  obtain  the  regionally  adjusted  or  lor  the  area 
being  studied,  a  data  set  used  to  evaluate  these  parameters  is  constructed  on  the 
basis  of  the  other  two  assumptions. 
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The  second  assumption  is  that,  at  short  epi central  distances,  equal  epicentral 
intensity  earthquakes  have  equal  levels  of  peak  acceleration,  regardless  of  regional 
source  variations.  Given  a  correlation  between  epicentral  intensity  and  accelera¬ 
tion  in  some  reference  region,  then  this  assumption  implies  that  the  relationship 
also  holds  in  any  other  region.  This  is  a  restricted  application  of  a  site  intensity - 
acceleration  correlation  in  that  only  the  more  impulsive  accelerations  are  assumed 
to  correlate  across  regional  boundaries.  In  this  way  the  regionally  dependent 
effects  of  distance  on  wavetrain  dispersion  are  either  eliminated  or  at  least  greatly 
reduced. 

The  use  of  epicentral  intensity,  as  opposed  to  magnitude,  as  the  basis  for 
inter-regional  comparison  of  near-field  accelerations  is  significant.  If  magnitude 
were  used  then  an  implied  assumption  would  be  made  that  earthquake  source  char¬ 
acteristics,  such  as  depth  of  occurrence,  are  regionally  invariant.  By  using  the 
surface  effects,  expressed  as  epicentral  intensity,  for  comparison  between  regions 
and  then  converting  intensity  to  magnitude  using  empirical  relations,  compensation 
for  these  source  variations  can  be  made. 

A  far-field  acceleration  estimate  can  be  made  by  assuming  that  the  peak 
acceleration  at  the  radius  of  the  felt  area  is  a  constant,  or  nearly  so.  The  basis 
of  the  third  assumption  is  that  the  minimum  acceleration  which  can  be  perceived 
by  a  human  is  a  function  more  of  the  human  anatomy  than  geology,  and  thus  should 
be  regionally  invariant.  The  radius  of  the  felt  area  as  a  function  of  epicentral 
intensity  can  be  obtained  from  historical,  noninstruniental  records.  Agai/i, 
empirical  relations  between  epicentral  intensity  and  magnitude  permit  the  conver¬ 
sion  of  this  information  to  a  magnitude  basis. 

The  question  then  arises  as  to  what  level  of  ac  celeration  is  appropriate  at  the 
limit  of  perceptibility.  Experiments  have  been  conducted  whic  h  have  found  that  the 
limit  of  perception  of  earthquake  type  motions  is  approximately  0.8  cm/sec •“  over 
the  frequency  range  of  2  to  0.4  Hz  (Ishimoto  and  Ootaka,  1022).  The  typical  per¬ 
ception  curve  which  was  found  is  shown  in  Figure  A2.  As  these  experiments  were 
conducted  under  laboratory  conditions  this  value  should  bo  considered  a  lower 
bound.  Field  experiments  conducted  at  the  same  time  found  a  threshold  value  for 
earthquake  perception  of  1.7  cm/s ee“  (Inouye,  1023).  Gultcnbcrg  and  Hichtcr 
(lOaG)  report  accelerations  over  the  range  of  1.6  to  2.0  cm/s  it'  for  Modified 
Merealli  intensities  of  I  and  II  in  California.  Based  on  the  definitions  of  the  inten¬ 
sity  scale,  the  limit  of  the'  felt  area  should  lie  between  these  two  intensities.  The 

consistency  of  these  values  suggest  that  a  reasonable  value  lor  use-  at  the  radius  of 

•> 

the'  fell  area  should  be-  between  1.7  and  4.0  rm/ser".  The'  higher  value  would  be 
more  applicable  in  a  region  of  low  seismic-  ae-Cvitv  where  the  association  of  slight 
vibrations  with  earthquakes  might  not  lie  made  as  readily  as  in  an  ac  tive  region. 

H2 


Figure  A2.  Minimum 
Perceivable  Accelerations 
lor  Periods  of  Less 
Than  2.  5  Seconds 


Period  (sec) 


(liven  these  assumptions  and  an  epieentrai  intcnsitv-acceleration  correlation 
for  the  reference  region,  a  data  set  ran  be  constructed  for  use  in  the  regression 
of  tin-  (y's  of  a  regionally  modified  attenuation  function.  The  data  set  would  consist 
of  a  near-field  acceleration  estimate  based  on  inter-regional  correlation  of  epi- 
rentral  intensity  and  a  constant  level  of  acceleration  at  the  radius  of  the  felt  area 
which  is  a  function  of  earthquake  magnitudes.  By  constraining  the  distance  factor, 

11  .  in  Lq.  (A6),  the  evaluation  of  the  constants  of  the  attenuation  function  can  be 
carried  out  itv  a  simple  least  squared  error  fitting  procedure. 

The  proposed  method  was  used  to  generate  aitenuation  functions  for  five  sub¬ 
divisions  of  tilt'  United  States.  These  regions  are  shown  in  Figure  A3  and  are  the 
same  as  those  used  bv  Brazee  (197  6)  to  develop  the  equations  relating  epieentrai 
intensity  and  magnitude  and  radius  of  felt  area  which  were  used  to  generate  t ho 
modified  attenuation  functions.  The  development  of  a  relation  for  California  w'as 
considered  a  partial  test  of  the  method.  This  application  could  only  test  the  ade¬ 
quacy  of  the  far-field  assumption  and  of  fitting  the  attenuation  function  to  a  limited 
data  base.  The  near-field  assumption  could  not  be  tested  as  the  available  epieentrai 
intensity -acceleration  correlations  are  based  on  California  data.  This,  however, 
is  not  considered  significant  as  all  other  procedures  for  obtaining  acceleration 
estimates  on  the  basis  of  site  intensity  make  use  of  a  more  liberal  form  of  ibis 
assum  ption. 

Before  proceeding  with  the  application,  two  points  should  lie  discussed.  The 
first  concerns  the  distance  at  which  the  near -field  estimate  of  acceleration  is 
made,  This  evaluation  c  an  be  made  either  at  some  fixed  distance  or  at  a  distance 
which  is  a  function  of  the  event  magnitude.  The  variable  distance  scheme  could 
result  in  such  large  distances  for  large  magnitude  events  that  dispersion  would 
become  an  important  factor  or,  for  small  magnitude  events,  such  small  distances 
that  anv  method  of  estimating  acceleration  would  he  invalid  because  of  very  near- 
field  effects.  Using  the  fixed  distance  requires  only  that  tilt'  distance  used  he 
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significantly  smaller  than  the  radius  of  tin*  felt  area  of  the  smallest  magnitude  used 
to  develop  the  data  set.  In  the  following  applications  of  the  method,  a  distance  of 
15  km  was  used. 


Figure  A3.  Regional  Divisions  l  sed  in  the 
Development  of  Derived  Acceleration  Attenua¬ 
tion  Functions 


Two  methods  are  also  available  for  the  evaluation  of  the  near-field  accelera¬ 
tions.  The  first  of  these  is  to  use  a  silt'  intensit v -acceleration  correlation  function 
of  the  form  of  Kq.  (A5)  or  to  convert  epicentral  intensity  to  equivalent  magnitude 
in  the  reference  region  and  to  use  a  reference  region  attenuation  function.  In 
general,  attenuation  functions  have  lower  standard  errors  than  available  site  in¬ 
tensity  relations,  hut  the  two  methods  are  probablv  equal  when  the  error  of  con¬ 
version  from  epicentral  intensity  to  magnitude'  is  included.  The  use  of  attenuation 
functions  was  selected  for  the  following  applications. 

California  -  The  basic  equation  for  evaluating  a  California  attenuation  function 
is  the  mnetion  defining  the  radius  of  the  felt  area,  K^,  in  terms  of  event  magnitude: 

In  U  0.  u.tti  .  o.  ;  15  M  .  t  A 7  > 

As  the  attenuation  lunction  used  to  evaluate  the  near-lield  accelerations  is  based 
on  data  from  California,  near-field  accelerations  can  he  determined  difeetlv  on 
tile  basis  ot  magnitude  w  ithout  the  intermediate  slen  of  epieent  ral  mfcusit  v  eomtiari 
son.  I  valuation  of  the  near-field  accelerations  were  made  using  two  differed 
acceleration  attenuation  functions: 


3.  10  ‘  0.8!'  ,M  -  1.  17  In  R 


and 


in  a  =  6.  16  *-  0.65  M  -  1.30  ln(K  *  25)  (A9) 

s 

both  of  which  were  developed  by  McGuire  (1978  and  1974,  respectively).  Both  2.0 
and  4.0  nn/sec'  were  used  as  the  acceleration  at  the  radius  of  the  felt  area.  The 
smallest  standard  error  of  fitting  were  found  using  Kq.  (A8)  and  4.  0  cm/sec". 

The  best  fit  equations,  based  on  the  standard  errors  of  fitting,  arc  given  by 

In  a  3.725  <  1.  026  M  -  1.  603  In  K  (A10) 

s 

and 

In  a  7.044  •  1.  155  M  -  2.  300  In  (R  •  25)  (All) 

s 

where  R  was  set  either  to  zero  or  25  km.  The  curves  generated  from  Kqs.  (A8) 
through  (All)  are  plotted  for  the  purposes  of  comparison  in  Figure  A4  evaluated  at 
5.0  M  and,  at  7.0  M  in  Figure  A5.  The  derived  curves  indicate  a  faster  fall  off  of 
peak  acceleration  than  either  of  the  empirical  functions.  This  becomes  most 
apparent  beyond  100  km.  If  the  far-field  assumption  is  correct,  this  difference  can 
probably  be  explained  bv  the  fact  that  the  empirical  equations  are  based  on  data 
predominately  from  shorter  epicentral  distances. 


Figure  A4.  Derived  (solid  curves)  and  Kmpiricul 
(dashed  curves)  Beak  Acceleration  Attenuation 
Functions  for  California  Fvaluated  at  5.0  M. 
Numbers  on  curves  refer  to  equation  numbers  givn 
in  text 
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Figure  A5.  Derived  (solid  curves)  and  Empirical 
(dashed  curves)  Peak  Acceleration  Attenuation 
Functions  for  California  Evaluated  at  7.0  M. 
Numbers  on  curves  refer  to  equation  numbers 
given  in  text 


In  Figure  A6,  all  four  equations  are  evaluated  at  6.  55  M  and  plotted  against 
sixty -one  rock  and  stiff  soil  peak  accelerations  recorded  during  the  Borrego 
Mountain  (6.5  M)  and  San  Fernando  (6.6  M)  Earthquakes  (Seed  et  al,  lfr6).  It  is 
apparent  the  curves  derived  bv  Ihe  method  proposed  in  this  report  do  a  better  job 
of  predicting  the  data  bevond  100  km.  The  mean  squared  residuals  of  the  natural 
logarithms  of  the  data  are  0.667,  0.587,  0.  lot,  and  0.220  |'or  Eqs.  (A8),  (A0), 

(A  10),  and  (All)  respective! v.  Inside  of  100  km  the  two  sets  of  equations  are 
essentially  equal  witli  mean  squared  residuals  between  0.  125  and  0.  105.  Outside 
of  this  distance  a  significant  change  occurs  with  mean  squared  residuals  of  0.  775, 
1.260,  0.  2)1-1,  and  0.  245,  respectively. 

Based  on  this  test,  it  would  appear  t lint  the  method  can  produce  a  reasonable 
attenuation  function  for  use  over  a  wide  range  of  distances.  In  addition,  the  succes 
of  the  I'ar-I'ield  assumption  would  suggest  that  it  could  be  used  as  an  additional 
restriction  in  the  development  of  empirical  functions  if  I'ar-I'ield  data  is  larking. 


Epicentra!  Distance  '  *rr 

Figure  A6.  Derived  (solid  curves)  and  Empirical 
(dashed  curves)  Peak  Acceleration  Attenuation 
Functions  for  California  Evaluated  at  6.  55  M,  and 
Plotted  Against  Rock  and  Stiff  Soil  Accelerations 
Observed  for  the  Borrego  Mountain  and  San  Fernando 
Earthquakes.  Numbers  on  curves  refer  to  equation 
numbers  given  in  text 

Other  Regions  -  Attenuation  functions  were  also  developed  for  four  other  re¬ 
gions  of  the  United  States.  The  relationships  used  in  the  development  of  these 
functions  and  the  resulting  functions  themselves  are  given  in  Table  Al.  The  de¬ 
rived  functions,  including  those  generated  for  California,  were  evaluated  for  a 
magnitude  of  6.  0  M  and  plotted  in  Figure  A",  for  those  with  Ro  set  to  0  and  in 
Figure  A8  for  those  with  R  i  set  to  25  km.  At  the  present  time  virtually  no  data 
exists  for  the  testing  of  these  functions.  The  functions  agree  with  Lhe  usually  hold 
concept  that  the  attenuation  throughout  the  western  United  States  is  fairly  uniform 
and  that  attenuation  is  much  lower  in  the  Central  United  States  with  the  Eastern 
United  States  being  somewhere  in  between. 

Finally,  a  comparison  of  the  accelerations  predicted  for  magnitude  4.0  M  and 
6.  0  M  earthquakes  in  the  Central  United  States  using  the  proposed  method  and  site 
intensity  based  methods  is  shown  in  Figure  AO.  The  site  intensity  methods  shown 
in  this  figure  use  the  intensit v -acceleration  and  intensity  and  distance-acceleration 
correlations  developed  bv  McGuire  (10771  in  conjunction  with  intensity  attenuation 
functions  developed  bv  Howell  and  Schultz  (1075).  While  it  is  impossible  al  the 
present  time  to  determine  which  of  the  equations  is  more  realistic,  il  is  apparent 
that  significant  differences  exist  between  then;. 


Acceleration  (cm/sec)  Acceleration  (cm/sec)  Acceleration  (cm /sec) 
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Figure  A?.  Derived 
Attenuation  Curves,  Using 
Rq  Set  to  Zero,  for  the 
Five  Subregions  of  the 
United  States  as  Given  in 
Table  Al.  All  curves 
evaluated  at  6.  0  M 


Figure  A8.  Derived 
Attenuation  Curves, 
Using  H0  Set  to  25  km, 
for  the  Five  Subregions 
of  the  United  States  as 
Given  in  Table  A  1. 

All  curves  evaluated 
at  6.  0  ]W 


Figure  A9.  Derived 
Attenuation  Curves  and 
Site  Intensity  [f(I)j  and 
Site  Intensity  and  Distances 
[  f(I,  R )]  Acceleration 
Attenuation  Curves  for  the 
Central  United  States. 

Kach  type  of  curve  is 
evaluated  at  4.  0  M  and 
6.  0  M 
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A4.  ERRORS 


It  is  not  possible  to  directly  evaluate  the  errors  associated  with  the  attenuation 
functions  derived  using  the  method  proposed  in  this  report.  One  scheme  which  does 
seem  reasonable,  however,  is  to  assume  that  the  distribution  of  accelerations  about 
the  derived  curves  is  similar  to  the  distribution  about  the  empirical  curves  for 
California.  A  standard  deviation  equal  to  those  evaluated  in  California,  or  perhaps- 
increased  to  be  conservative,  can  be  assigned  to  these  functions.  A  typical  standard 
deviation,  based  on  the  values  listed  by  McGuire  (1976)  is  about  0.707  for  the 
natural  logarithm  of  the  peak  acceleration  and  might  be  appropriate  for  use  with 
the  regionally  modified  attenuation  functions. 


A5.  CONCLUSIONS 

Although  there  are  many  significant  problems  associated  with  the  use  of  strong 
ground  motion  attenuation  functions,  they  do  provide  an  important  means  for  con¬ 
verting  knowledge  of  seismic  activity  into  estimates  of  engineering  seismic  risk. 

In  this  report  a  new  method  for  the  regional  modification  of  peak  acceleration 
attenuation  functions  has  been  presented  and  the  method  applied  to  derive  attenua¬ 
tion  functions  for  five  regions  in  the  United  States.  The  method  relies  on  two 
primary  assumptions.  First,  in  the  near-field,  the  level  of  acceleration  causing 
a  given  epicentral  intensity  is  regionally  invariant.  Second,  that  the  level  of 
acceleration  at  the  radius  of  the  felt  area  is  a  constant.  Available  data  tends  to 
support  the  second  assumption  while  the  first  is  a  restricted  application  of  the 
assumption  used  in  previously  suggested  methods. 

The  application  of  this  technique  to  the  development  of  an  attenuation  function 
for  California  indicates  the  feasibility  of  the  approach.  The  derived  functions 
predict  peak  accelerations  for  the  Borrego  Mountain  and  San  Fernando  Earthquakes 
at  least  as  well  as  the  two  empirical  equations  to  which  they  were  compared.  Be¬ 
yond  100  km,  the  derived  functions  were  found  to  be  significantly  more  accurate 
than  the  empirical  functions. 

Peak  acceleration  attenuation  functions  were  also  generated  for  four  other 
regions  of  the  United  States.  While  the  errors  associated  with  these  functions  can 
not  be  directly  calculated,  it  is  suggested  that  the  errors  might  be  similar  to  those 
found  for  empirical  functions.  One  advantage  oT  the  proposed  method  is  that,  as 
observational  data  becomes  available  in  the  region  for  which  the  modification  is 
conducted,  it  can  be  incorporated  into  the  data  base  used  lor  regression  of  the 
attenuation  function  parameters,  thus  providing  an  improving  estimate  of  the 
regional  attenuation  function  until  the  point  is  reached  where  the  observational 
data  alone  is  sufficient. 
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